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FOREWORD 


The  calorimetry  of  any  electrochemical  cell  involves  two  type  of  activities:  data 
collection  and  data  evaluation.  The  required  data  are  the  cell  potential-time 
and  cell  temperature-time  series.  The  evaluation  is  based  on  conservation  laws 
subject  to  constraints  dictated  by  cell  design  and  the  adapted  experimental 
procedure. 

Volume  2  of  this  report  deals  with  the  modeling  and  simulation  of  the  Dewar-type 
calorimeter.  It  was  written  by  Professor  Fleischmann  to  provide  an  authoritative 
discussion  of  the  calorimetry  of  electrochemical  cells.  The  emphasis  is  on  the 
interpretation  of  data  and  the  accuracy  of  the  determination  of  the  excess  en¬ 
thalpy  generation  via  the  appropriate  selection  of  heat  transfer  coefBcients.  The 
discussion  of  the  calorimetry  of  the  Dewar-type  cells  is  presented  in  the  form  of 
technical  report  for  a  number  of  reasons,  among  them:  (i)  its  length  would  likely 
prohibit  publication  in  topical  journals,  (ii)  to  clarify  misunderstandings  regard¬ 
ing  the  principles  of  calorimetry  as  applied  to  electrochemical  cell  in  general  and 
to  the  cell  employed  by  Fleischmann  and  his  collaborators,  in  particular. 


S.  Szpak  and  P.A.  Mosier-Boss,  eds. 
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INTRODUCTION 

Apart  from  some  fragmentary  investigations,  primarily  related  to  the  study  of 
the  self“discharge  of  batteries,  there  exists  no  well  defined  set  of  studies  in  the 
field  of  the  electrochemical  calorimetry.  We  note  that  such  studies  would  allow 
the  investigation  of  the  thermal  behavior  of  a  wide  range  of  reactions,  especially 
irreversible  processes.  Thus,  the  establishment  of  an  accurate  model  of  an  ex¬ 
periment  is  very  important.  However,  as  this  aspect  is  not  generally  understood, 
we  felt  it  necessary  to  produce  this  document. 

In  spite  of  its  length,  this  volume  only  covers  the  analysis  of  a  data  set  generated 
by  calculation  and  one  measurement  cycle  for  a  “blank  experiment,”  We  believe 
that  it  is  very  important  to  produce  a  detailed  analysis  and  account  (as  far 
as  is  possible  at  this  stage)  of  the  methodology  which  we  adopted.  This  is 
especially  important  in  view  of  the  misleading  comments  which  have  been  made 
about  the  calorimetry  of  the  Pd/D  system.  Taken  at  face  value,  one  must 
believe  that  the  workers  concerned  do  not  understand  the  difference  between 
differential  and  integral  coeflicients,  the  disadvantages  of  differentiating  “noisy” 
data  as  compared  to  integrating  such  data,  the  differences  between  the  precision 
and  accuracy  of  data  evaluations,  the  recognition  of  “negative”  and  “positive 
feedback,”  the  analysis  of  cooling  curves,  etc.  They  do  not  understand  relaxation 
nor  recognize  the  presence  of  strange  attractors  and  the  way  in  which  the  effects 
of  such  complications  can  be  circumvented.^ 

It  is  relevant  here  to  reflect  on  the  precision  and  accuracy  of  the  experiments. 
Of  course,  if  the  precision  is  high,  then  there  will  be  no  difficulty  in  interpreting 
changes  in  the  rates  of  excess  enthalpy  generation  as  small  as  1  mW  at  the  10  a 

^Of  course,  it  is  possible  that  the  researchers  concerned  do  not  understand  any  of  these 
matters,  but  what  is  so  remarkable  is  that  they  have  failed  to  understand  these  topics  even 
when  they  have  been  described  to  them. 
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level.^  Of  course,  the  question  of  the  magnitude  of  the  errors  raises  three  further 
important  questions:  (i)  what  error  limits  are  required  so  as  to  be  able  to  detect 
excess  enthalpy  generation  at  an  adequate  level  of  statistical  significance?  (ii) 
what  is  the  difference  (if  any)  between  the  experiments  carried  out  with  ICARUS 
systems  and  ICARUS  look-alikes  and  with  other  types  of  calorimetry?  (iii)  how 
can  one  assess  the  error  limits  of  a  given  piece  of  instrumentation? 


The  answer  is  that  one  simply  stops  the  development  of  the  methodology  when 
one  is  able  to  make  an  adequate  set  of  measurements.  We  note  here  that  this 
particular  specification  is  itself  dependent  on  the  physical  size  of  the  systems 
being  investigated  as  well  as  the  chosen  operating  conditions.  In  our  particular 
investigation  the  limit  was  certainly  reached  when  the  errors  had  been  reduced 
to  the  0.01%  level.  Naturally,  the  first  question  impacts  on  the  second  and  we 
note  that  it  is  the  use  of  less  precise  and  accurate  calorimetric  methods  which  has 
bedeviled  so  much  of  the  research  in  this  field.  The  reason  is  that  with  the  use 
of  less  precise/accurate  methods,  it  becomes  impossible  to  monitor  the  build-up 
of  excess  enthalpy  generation.  This  then  brings  us  to  the  third  question  and  the 
answer  to  this  is  exactly  with  the  methods  outlined  in  this  document,  at  least 
as  far  as  isoperibolic  calorimetry  is  concerned  (although  it  is  not  very  difficult 
to  specify  improvements  in  those  methods!).^  It  is  relevant  that  although  errors 
had  undoubtedly  been  made  in  setting  up  these  experiments,  the  detailed  data 
analyses  had  also  shown  the  way  in  which  such  errors  could  be  allowed  for.^ 


To  reiterate,  we  considered  it  necessary  to  produce  this  document  for  the  follow¬ 
ing  reasons:  Firstly,  it  is  always  essential  to  determine  the  Instrument  Function 
(or  of  a  parameter  or  sets  of  parameters  which  define  the  Instrument  Function) 
and  to  validate  the  methods  of  data  analysis.  Such  validation  is  best  done  us¬ 
ing  simulated/calculated  data.  Secondly,  one  needs  to  see  the  extent  to  which 
“blank”  experiments  conform  to  expectations.  Thirdly,  one  needs  to  investigate 
the  ways  in  which  methods  of  data  analysis  may  fail. 

^However,  the  high  precision  of  the  instrumentation  (relative  errors  below  0.01%)  has  been 
converted  into  a  10%  error  by  the  group  at  NHE.  It  is  hard  to  see  how  anybody  can  make  such 
an  assertion  while  still  keeping  a  straight  face.  If  the  errors  were  as  high  as  this,  then  it  would 
be  impossible  to  say  anything  sensible  about  calorimetry  -  for  that  matter,  it  would  remove 
one  of  the  main  planks  of  scientific  methodology. 

^The  answer  to  this  question  brings  us  to  very  interesting  further  lines  of  enquiry  which 
can  be  summarized  by  the  question:  “why  is  it  that  NHE  have  never  made  any  sets  of  raw 
data  for  blank  experiments  available  for  further  analysis?”  If  one  considers  this  question  in  a 
naive  way,  then  one  would  say  that  there  can  hardly  be  any  reason  for  not  releasing  data  sets 
which  do  not  show  any  generation  of  excess  enthalpy! 

"^Instead  of  seeking  to  establish  the  correct  way(s)  of  calibrating  the  systems,  the  group 
at  NHE  used  the  procedure  leading  to  (A:'^°)362,  probably  coupled  to  timing  errors  in  the 
calibration  pulse  which  they  did  not  allow  for.  Needless  to  say,  this  produced  nonsensical 
results  which  they  used  as  a  justification  for  substituting  an  invalid  method  of  data  analysis. 
Moreover,  this  invalid  method  of  data  analysis  was  applied  to  just  two  experiments,  regarded 
as  being  typical,  although  the  fact  that  there  were  malfunctions  in  these  experiments  has  also 
been  pointed  out. 
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SYMBOLS  USED 

C  -  heat  capacitance.  [J(gMole)"^K"^] 

Ec[t)  -  cell  voltage  at  time  t.  [V] 

Eth,b  ~  thermoneutral  potential  at  bath  temperature.  [V] 

F  -  Faraday  constant,  [coulombs (gMole)^] 

AH  -  rate  of  enthalpy  input.  [W] 

AH^^  -  rate  of  evaporative  cooling.  [W] 

AHnet{i)  -  rate  of  net  enthalpy  input  at  time  t.  [W] 
k  -  heat  transfer  coefficient.  [WK“^] 

L  -  latent  heat  of  evaporation.  [J(gMole)“^] 

M  “  number  of  moles  of  D2O  at  t  =  0. 

P  -  vapor  pressure  at  the  cell  temperature,  [bar] 

P*  -  atmospheric  pressure,  [bar] 

Qf{t)  “  rate  of  generation  of  excess  enthalpy  in  the  cell  at  time  t.  [W] 
t  -  time,  [s] 

a  -  defined  in  Eq.  (27). 

7  -  related  to  time  dependence  of  the  (k)  coefficient  defined  in  Eq.  (20). 
A  -  defined  in  Eq.  (35). 

A  -  defined  in  Eq.  (36). 

T  -  time,  [s] 

0  -  bath  temperature.  [K] 

-  temperature  difference  between  the  cell  and  the  water  bath.  [K] 
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SECTION  1:  THE  EVOLUTION  OF  THE  ICARUS  DATA  EVALU¬ 
ATION  STRATEGIES. 


We  have  in  the  past  used  a  variety  of  strategies  to  evaluate  the  experimentally 
determined  temperature-time  and  cell  potential-time  series.  These  strategies 
can  be  described  (at  least  in  part)  by  the  heat  transfer  coefficients  which  govern 
the  “Black  Box”  representing  the  calorimeters.  We  will  confine  attention  here 
to  the  particular  forms  of  the  heat  transfer  coefficients  which  were  important  for 
the  evolution  of  the  data  evaluation  strategies  used  with  the  ICARUS  Systems. 
There  is  a  considerable  amount  of  material  which  needs  to  be  considered  even  if 
we  pose  this  restriction. 

The  specification  of  the  ICARUS  Data  Evaluation  Strategies  [1]  (as  modified  in 
part  in  [2])  was  based  mainly  on  the  analysis  of  temperature-time  and  input 
enthalpy-time  series  generated  by  simulations  based  on  the  differential  equa¬ 
tions  representing  the  models  of  the  calorimeters  as  well  as  on  the  evaluation 
of  suitable  “blank  experiments”  (principally  using  Pt  cathodes  polarized  in  0.1 
M  Li0D/D20),  compare  [3-6]).  For  the  first  of  these  sets  of  tests,  the  validity 
(or  otherwise)  of  the  methods  used  was  judged  by  the  recovery  of  the  parame¬ 
ters  used  in  the  simulations  at  adequately  high  levels  of  statistical  significance. 
These  tests  have  now  been  reconstructed  and  are  described  in  Section  5  where 
it  will  be  seen  that  they  were  sufficient  for  the  specification  of  the  ICARUS  data 
processing  strategies  (as  well  as  for  the  specification  of  some  of  the  shortcomings 
of  the  various  methodologies). 


The  validity  (or  otherwise)  of  the  methods  defined  in  the  first  set  of  tests  was 
then  further  evaluated  by  using  data  for  “blank  experiments.”  This  validity  was 
then  assessed  principally  by  investigating  the  degree  of  conformity  with  the  pre¬ 
dictions  based  on  the  first  set  of  tests  as  well  as  by  the  statistical  significance  of 
the  methodologies;  by  determining  the  degree  of  correspondence  of  the  “true” 
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and  “lower  bound”  heat  transfer  coeflBcients  and,  related  to  this,  the  observation 
of  a  “zero”  rate  of  excess  enthalpy  generation  [4];  by  checking  the  relationship  be¬ 
tween  the  various  forms  of  the  heat  transfer  coefficient  [3^];  and  by  determining 
whether  the  experimentally  observed  relaxations  of  the  temperature-time  series 
conform  to  predictions  based  on  simplified  models  of  the  calorimeters  [4].  One 
such  set  of  tests  for  a  “blank”  experiment  carried  out  using  an  ICARUS-2  sys¬ 
tem  is  described  in  Section  6  (but  excluding  the  test  of  the  relaxation  behavior). 
The  conclusions  drawn  from  tests  such  as  those  described  in  Sections  5  and  6 
are  summarized  in  Section  7,  where  they  are  compared  with  the  specifications 
of  the  experimental  protocols  and  data  evaluation  strategies  contained  in  the 
Handbook  for  the  ICARUS-1  Systems  [1]. 

It  is  perhaps  not  surprising  that  most  of  this  material  has  never  been  pub¬ 
lished.  We  could  not  imagine  that  any  paper  written  on  this  topic  would  ever 
be  accepted  by  a  scientific  journal  and  believed  that  it  would  be  sufficient  for 
us  to  specify  the  protocols  and  data  evaluation  strategies  to  be  used  with  the 
ICARUS  Systems  [1].  It  is  the  deviation  of  the  NHE  group  from  the  recom¬ 
mended  protocols  and  their  use  of  inaccurate  and/or  invalid  methods  of  data 
evaluation  [7]  which  makes  it  necessary  to  reconsider  the  background  material 
for  the  ICARUS- 1  systems.  Unfortunately,  such  a  reconsideration  makes  it  nec¬ 
essary  to  give  a  more  closely  defined  description  for  some  of  the  heat  transfer 
coefficients  than  has  hitherto  been  used.  This  is  contained  in  Section  2.  We  will 
use  the  designation  where  i  —  1,2,3  denotes  “differential”,  “backward 

integration”  and  “forward  integration”  respectively;  j  is  defined  in  Section  2  and 
Z  —  1,2  denotes  “lower  bound”  and  “true”  respectively.  We  had  hoped  to  cir¬ 
cumvent  the  need  for  such  an  extended  description  so  as  to  avoid  overburdening 
our  accounts  with  redundant  symbolism.  Evidently,  we  were  mistaken  with  our 
descriptions. 

Section  4  describes  the  method  used  in  the  simulations  of  the  data  evaluated  in 
Section  5.  It  will  be  seen  that  the  first  set  of  tests  were  incomplete  although 
they  were  sufficient  for  the  specification  of  the  second  set  of  tests  (Section  6) 
and,  in  turn,  for  the  specification  of  the  ICARUS  protocols  and  data  processing 
strategies  [1], 
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SECTION  2:  DEFINITION  OF  THE  HEAT  TRANSFER  COEFFI¬ 
CIENTS. 


It  is  convenient  to  consider  the  form  of  the  cell  temperature-time  series  generated 
by  the  simulation  described  in  Section  3,  Fig.  1.  As  will  be  shown  in  that  Section, 
the  data  have  been  generated  using  Eq.  (31)  in  which  the  heat  transfer  from 
the  cell  is  described  as  being  pseudo  -  conductive  while  the  effects  of  “negative 
feedback”  are  taken  into  account  by  writing  the  input  enthalpy  as 

input  enthalpy  =  (E  —  AA^)I,  (1) 

where 

E  =  Ee{t  =  0)-Eth,t  (2) 

and 

A0{t)  =  d{t)-et.  (3) 

However,  in  order  to  prepare  the  way  for  the  consideration  of  the  “blank  ex¬ 
periment”,  Section  5  as  well  as  the  later  sections,  the  heat  transfer  coefficients 
are  defined  by  describing  the  heat  transfer  as  being  pseudo-radiative,  Eq.  (21), 
giving  the  “temperature  function” 

Me)^[9i  +  A9it)]^-0t-  (4) 

For  the  case  of  the  description  of  heat  transfer  being  pseudo-conductive,  we  need 
to  replace  the  temperature  function  by  A0, 

h{e)  =  Ae,  (5) 


We  can  then  define  a  “lower  bound  heat  transfer  coefficient”  (i.e.,  a  coefficient 
which  assumes  that  the  rate  of  excess  enthalpy  generation  is  zero)  for  any  part 
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(6) 


of  the  measurement  cycle,  Fig.  1, 


[Ecit)  -  Eth,t]I  -  -  CpMidAO/dt) 

fm 


In  our  early  work  we  evaluated  this  coefficient  just  before  the  start  of  the  cal¬ 
ibration  pulse  and  designated  this  particular  value  as  (A;^)4  (see  [8]),  a  matter 
which  was  evidently  not  understood  [9]  although  the  calculation  scheme  was  set 
out  in  Appendix  4  of  [8].  The  special  value  of  (fc^)ii  just  before  the  end  of  the 
calibration  pulse  had  originally  been  designated  as  (fc^)i.  (Together  with  (A:^)2 
described  below,  these  were  the  first  two  coefficients  which  we  used  in  the  data 
analysis,  hence  their  designation  with  the  suffixes  1  and  2).  It  was  our  investiga¬ 
tion  of  the  “Harwell  Data  Sets”  [10]  which  convinced  us  that  the  “lower  bound 
heat  transfer  coefficient”  is  more  useful  than  just  the  two  special  values  (/r^)i 
and  (fc/j)4  leading  to  (A;^)ii  as  derived  in  the  “ICARUS  (fc/H)ii -spreadsheets,” 
see  Section  6. 

Having  obtained  (A:^)ii,  we  frequently  wish  to  establish  the  11-point  averages 
(fci^)ii  so  as  to  decrease  the  “noise.”  This  gives  us  ca  26  values  for  measurement 
cycles  lasting  1  day  or,  better  ca  52  values  for  the  recommended  2-day  cycles. 
In  turn  it  is  useful  to  evaluate  the  6-point  averages  of  (fci^)ii  which  we  have 

designated  as  It  is  not  useful  to  extend  this  averaging  beyond  6  points 

because  any  such  extension  makes  the  systematic  errors  (due  to  the  residual 
decrease  of  with  time)  larger  than  the  random  errors,  that  is,  if  the 

systems  are  behaving  sensibly. 


It  will  become  apparent  that  we  need  accurate  values  of  CpM  to  make  (A:^)ii 
generally  useful  but,  if  we  exclude  regions  where  the  temperature  is  varying 
rapidly  with  time,  then  “guesstimates”  of  CpM  are  quite  adequate.  We  note 
that  if  we  rearrange  Eq.(6)  to  the  straight  line  form 

y  —  mx  -h  c,  (7) 


i.e., 


[Ea{t)  -  Eth,bV  ^  _  CpMjdAe/dt) 


fiiO) 


h{0) 


(8) 


then  approximate  values  of  CpM  can  be  obtained  from  the  slopes  of  the  plots  de¬ 
rived  for  regions  where  the  temperature  is  varying  relatively  rapidly  with  time. 
We  can  distinguish  four  such  plots  which  we  have  designated  by  the  relevant 
derived  heat  transfer  coefficients:  (A:^^)i5i,  (^i?^)i7i  and  (fc^^)i8i  ac¬ 

cording  to  whether  the  fitting  of  Eq.  (8)  is  carried  at  at  times  somewhat  above 
the  origin,  at  times  somewhat  above  ti  (the  time  of  application  of  the  calibra¬ 
tion  pulse) ,  at  times  somewhat  above  t2  (the  time  of  cessation  of  the  calibration 
pulse)  or  by  the  combination  of  the  last  two  time  regions.  Fig.  1.  However,  we 
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note  that  there  is  a  measure  of  ambiguity  about  the  interpretation  of  the  values 
of(A:^)i,j,i  derived,  which  will  be  discussed  in  Sections  5  and  6. 

We  note  here  that  separate  investigations  showed  that  (dA0(i)/dt)  is  best  esti¬ 
mated  by  using  the  second  order  central  differences  (i.e.,  the  chords  of  the  curves) 
when  using  “real”  data  (i.e.,  experimental  rather  than  simulated  data).  More 
accurate  values  could  be  derived  in  principle  by  using  higher  order  differences. 
However,  in  practice,  the  repeated  differentiation  of  the  experimental  data  leads 
to  an  increase  in  “noise”  if  we  use  differences  higher  than  the  second  order.  This 
use  of  the  central  difference  is  of  some  importance  when  we  consider  the  integra¬ 
tion  processes  required  for  the  derivation  of  the  heat  transfer  coefficients  based 
on  the  forward  integration  of  the  experimental  data  (see  below) . 

We  note  also  that  objections  have  often  been  raised  to  the  procedures  which  we 
have  adopted  based  on  the  fact  that  we  have  not  “binned  the  data,”  i.e.,  we  have 
not  signal  averaged  before  the  data  analysis.  However,  “binning  of  the  data” 
must  always  be  approached  with  great  caution:  one  should  only  “bin  data”  or 
“bin  coefficients”  if  these  data  or  coefficients  can  be  expected  to  be  constant 
over  the  averaging  interval.  This  is  not  the  case  for  (fc^)ii  unless  the  effects  of 
the  term  CpM{dM{t)ldt)  have  been  taken  into  account.  Once  this  is  done,  we 

can,  of  course,  bin  the  data  as  we  have  done  in  deriving  and  (^^2)11* 


In  the  case  of  the  interpretation  of  data  derived  with  calorimeters  relying  on 
radiative  cooling,  the  position  is  further  complicated  by  the  fact  that  the  differ¬ 
ential  equation  representing  the  calorimeters,  Eq.  (20),  is  both  nonlinear  and 
inhomogeneous.  It  does  not  follow  therefore  that  coefficients  derived  by  averag¬ 
ing  the  data  are  the  same  as  averages  of  the  coefficients  derived  by  using  the  raw 
data.  However,  we  did  confirm  in  1992  (when  this  whole  saga  was  first  investi¬ 
gated)  that  we  do  in  fact  obtain  an  equivalence,  provided  we  restrict  attention 
to  regions  where  {dA0{t)/dt)  is  adequately  small.  We  concluded  that  in  such 
regions  the  differential  Eq,  (20)  could  be  sufficiently  linearized  and  that  second 
order  small  differences  were  sufficiently  small  to  allow  such  averaging.  However, 
as  the  values  of  (/c/^)ii  obtained  in  this  way  were  identical  to  those  obtained 
following  the  procedure  outlined  above,  there  was  evidently  no  justification  in 
pursuing  the  matter  further  (nor  to  complicate  the  instrumentation!). 


We  will  refer  to  the  averaging  procedures  further  below  when  discussing  the 
heat  transfer  coefficients  based  on  the  forward  and  backward  integration  proce¬ 
dures,  It  is  next  necessary  to  evaluate  the  “true  heat  transfer  coefficients.”  The 
value  {k^fi)2  near  the  end  of  the  calibration  period  is  obtained  by  including  the 
calibration  pulse,  AQ: 
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(^r)2 


AQ  +  [£^c(^^2! ^2)  ~  Ec{AOi ,  ti]I  —  AHev {A92, H)  +  AH^viAdi ,  ti) 

MO) 


C,M[idAe/dt)A0,,,,  -  d{A9/dt)Ae,M 

MS) 

where  we  now  have 


MS)  =  [^6  +  (A02,tJ]'’  -  [06  +  (A0i.tj^ 


(9) 

(10) 


In  order  to  carry  out  this  evaluation,  it  is  useful  to  construct  A4-  or  A3-  sized 
plots  (European  notation);  see  Sections  5  and  6,  and  then  to  obtain  appropriate 
averages  using  a  transparent  ruler.  This  type  of  analysis  used  to  be  a  generally 
accepted  approach  but  then  fell  into  disrepute.  However,  it  is  now  again  accepted 
as  giving  so-called  “robust  estimates.” 


We  note  that  the  errors  in  (/i:^)2  are  measures  of  the  accuracy  of  the  “true 
heat  transfer  coefRcient”  as  the  estimate  is  made  in  terms  of  the  known  Joule 
enthalpy  input  to  the  calibration  heater.  Errors  in  (A:^)i  or  (/:ii)ii  are  measures 
of  the  precision  of  the  “lower  bound  heat  transfer  coefficients”  as  there  is  no 
independent  calibration  and  there  may  be  excess  enthalpy  generation  in  the 
system.  It  is  important  that  (fc^)ii  and  {k'^)2  are  the  least  precise  and  least 
accurate  coefficients,  which  can  be  obtained  from  the  raw  data.^ 


We  have  always  insisted  that  the  construction  and  evaluation  of  plots  of  the 
raw  data  is  an  essential  prerequisite  of  the  more  elaborate  data  evaluation  pro¬ 
cedures.  For  one  thing,  it  shows  whether  the  “noise  levels”  in  the  experiments 
were  sufficiently  low  to  justify  more  detailed  evaluations  and  also  points  to  mal¬ 
functions  in  the  operation  of  the  experiments.  It  is  very  important  therefore 
to  establish  whether  the  group  at  NHE  ever  followed  this  particular  instruction 
and,  if  they  did,  what  conclusions  they  may  have  drawn  from  any  such  plots. 


Having  obtained  the  “true  heat  transfer  coefficient”  at  a  single  point  (usually 
near  the  end  of  the  calibration  pulse),  it  is  important  to  ask:  “what  is  the 
true  heat  transfer  coefficient  (fe^)i2  at  any  other  time?”  We  can  make  such  an 
evaluation  within  the  duration  ti  <  t  <  of  the  calibration  pulse  simply  by 
using  Eq.  (9),  giving  (A;^)i2  rather  than  (A:^)2-  Note  also  that  Eq.  (9)  can  be 
rearranged  to  the  straight  line  form 

®Any  statements  that  the  errors  are  larger  than  this  (as  has  been  made,  for  example,  in 
the  paper  from  the  group  at  NHE  [7])  simply  show  that  mistakes  have  been  made  in  the  data 
analysis  procedures  and/or  the  execution  of  the  experiments. 
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AQ  +  [Ec{Ad2,t)  -  Ec{Aei,t)]I  -  AHev{Ae2,t)  +  AHev{A9i,t) 

/2W 


+  {k'^)lG2y 


(11) 


_  CpM[{dAe/dt)Ae2.,  -  {dA9/dt)Ae^,A 

12(0) 

which  is  applicable  to  times  close  to  and  above  ti .  It  is  evident,  therefore,  that 
such  data  derived  from  the  experiments  can  also  be  used  to  obtain  estimates  of 
CpM,  but  the  accuracy  of  such  values  is  inevitably  much  lower  than  of  those 
obtained  by  the  application  of  the  corresponding  expression  for  the  “lower  bound 
heat  transfer  coefficient,”  (fc^^)i6i5  Eq.  (8).  Nevertheless,  Eq.  (11)  is  useful 
because  it  allows  the  removal  of  the  effects  of  the  water  equivalent  on  the  “true 
heat  transfer  coefficient,”  (/c^°)i62,  simply  by  extrapolating  to  zero  value  of  the 
abscissa. 


In  the  regions  in  which  there  is  no  application  of  the  heater  pulse,  i.e.,  for 
0  <  t  <  h  and  t2  <  t  <  T,  the  “true  heat  transfer  coefficient”  can  only  be 
obtained  from  the  “heating”  and  “cooling  curves,”  i.e.,  the  “driving  force”  is  the 
change  in  the  enthalpy  content  of  the  calorimeters  rather  than  AQ.  It  is  now 
sensible  to  cast  Eq.  (9)  in  the  form 


Cj,M[{dAe/dt)Ae,,t  -  {dAe/dt)Ae,, 

hie) 


-  -(4°) 


152 


[E,(Ag2,0  -  E,{Mut)]I  -  +  AHef  (Agl,f) 

f2{e) 


(12) 


If  the  system  is  functioning  correctly,  then  it  will  be  found  that  the  L.H.S.  of  Eq. 
(12)  is  essentially  constant  (although  this  constancy  can  only  be  probed  over  a 
short  time  range).  The  second  term  on  the  R.H.S.  of  Eq.  (12)  will  be  much 
smaller  than  the  term  on  the  L.H.S. ,  i.e.,  it  is  in  the  nature  of  a  correction  term 
to  give  the  “point-by-point”  values  of  (/!:^^)i52  or  (fc^^)i72*  It  will  be  evident 
that  the  accuracy  of  these  versions  of  the  “true  heat  transfer  coefficient”  is 
limited  by  the  accuracy  of  our  estimates  of  CpM.  This  particular  part  of  the 
methodology  is  therefore  only  useful  to  serve  as  a  check  on  the  operation  of  the 
cells  and  methods  of  data  evaluation.®  The  importance  of  Eqs.  (11)  and  (12) 
lies  partly  in  the  fact  that  these  equations  lead  to  the  interpretation  of  what 
we  have  classified  as  Cases  1  and  2  of  the  phenomenon  of  “Heat  after  Death”  [3] 
(and  to  other  related  versions  of  these  Cases). 

The  assumption  underlying  this  part  of  the  account  is  that  we  can  only  determine 
(^k)i2  within  the  duration  of  application  of  the  calibration  pulse, 

is  not  possible  to  combine  the  data  in  the  regions  just  above  t\  and  t2  to  give  a  simple 
equation  leading  to  at  least  we  should  say  that  we  have  not  been  able  to  specify  a 

simple  data  processing  strategy! 
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t\  <  t  <  h,  Fig  1,  and  at  lower  accuracy,  and  {k'^)n2,  in  regions 

adjacent  to  the  origin  and  for  times  adjacent  and  above  h-  However,  this  con¬ 
clusion  is  incorrect.  We  need  to  make  the  additional  assumption  that  the  rate 
of  any  excess  enthalpy  generation  is  constant  during  any  particular  calibration 
period  in  order  to  determine  (fcfi)i2-  This  means  that  we  can  only  obtain  a  single 
value  of  this  heat  transfer  coefficient  per  calibration  period  and,  consequently, 
a  single  value  of  (fcfl)i2  -  (fcfl)ii-  Two  important  points  follow  from  this  con¬ 
clusion.  In  the  first  place,  the  precision  of  must  be  very  nearly  equal  to 

the  precision  of  {k'fi)i\  (this  is  discussed  further  in  Section  6).  Secondly,  if  we 
extend  the  assumption  that  the  rate  of  excess  enthalpy  generation  is  constant 
during  the  period  ti  <  t  <  t2  to  saying  that  it  is  constant  during  the  period 
0  <t  <T,  then  it  is  immediately  possible  to  derive  (A:/j)i2  over  the  whole  of  the 
measurement  cycle.  The  values  obtained  can  be  compared  to  the  special  values 
of  {k'^)i52  and  (A;^°)i72  in  the  relevant  time  regions. 

Having  obtained  the  local  (differential)  values  of  the  “lower  bound”  and  “true 
heat  transfer  coefficient,”  we  naturally  cast  around  for  methods  which  would 
increase  their  precision  and  accuracy.  The  reason  for  the  limited  precision  and 
accuracy  is  mainly  due  to  the  need  to  differentiate  the  noisy  experimental  data 
sets.  In  our  early  work,  we  overcame  this  particular  difficulty  by  using  {kji)4 
and  (fcfl)2  as  starting  values  for  the  nonlinear  regression  procedure  leading  to 
the  heat  transfer  coefficient  (/cijjs.  Here,  we  fitted  the  numerical  integrals  of  the 
differential  equation  governing  the  behavior  of  the  calorimeters  to  the  data  sets 
for  complete  measurements  cycles.  The  relative  errors  in  (fcjjjs  which  we  could 
achieve  in  this  way  were  below  0.1%,  even  when  using  the  unsatisfactory  early 
version  of  our  calorimeters  (i.e.,  those  not  silvered  in  their  top  portions). 

The  use  of  nonlinear  regression  procedure  had  the  further  distinct  advantage  that 
we  could  adjust  the  integration  interval  in  regions  where  the  temperature  was 
varying  rapidly  with  time  so  as  to  achieve  the  required  accuracy  in  the  integrals. 
This  is  not  possible  for  the  methods  which  we  will  discuss  below  (and  which 
were  part  of  the  ICARUS  data  processing  strategy)  because  the  intervals  for  the 
data  acquisition  were  fixed.  As  a  matter  of  fact,  the  interval  300  s  was  chosen 
because  such  an  interval  does  not  degrade  the  evaluation  of  any  of  the  series 

{k'ji)21,  {k'ji)22,  ik'i^)251y  {k'^)2o2,  ik'^)26\ ,  {k'^)262,  (fcfl°)271>  and  (fc^°)272- 

However,  it  does  degrade  the  evaluation  of  ik'n)3i,  (A:^°)35ii  (A:'’°«)36i,  and 
{k'^)37i  to  some  extent,  and  leads  to  a  marked  degradation  of  {k'n)32,  {k'^)332, 
{k'^)362,  and  {k'^)372-  The  fact  that  the  data  acquisition  interval  was  too  long 
for  straightforward  estimations  of  the  (A:jj)3,j,j  series  of  heat  transfer  coefficients 
was  already  pointed  out  to  NHE  in  the  ICARUS-1  Handbook  [1].  These  matters 
will  be  considered  further  in  Sections  5  and  6. 

We  have  pointed  out  on  other  occasions  that  the  reason  we  opted  for  using 
nonlinear  regression  fitting  in  our  early  work  was  because  the  pressure  of  events 
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did  not  allow  us  to  go  through  the  logical  sequence  of  using  linear  regression, 
multilinear  regression  and  nonlinear  regression  (in  fact,  we  had  to  opt  for  a 
“catch-all”  methodology).  However,  as  we  could  not  make  nonlinear  regression 
“user  friendly”  with  the  computing  power  then  available  to  us,  so  in  1991-92,  we 
investigated  the  application  of  linear  regression,  which  became  part  and  parcel 
of  the  ICARUS- 1  methodology.*^ 


Attention  has  been  drawn  to  some  of  these  details  because  it  would  certainly  be 
possible  to  reimplement  parts  of  these  projects  provided  one  could  gain  access 
to  the  data  sets.  If  one  wishes  to  avoid  the  numerical  differentiation  of  the 
experimental  data  sets,  then  one  can  rely  instead  on  the  numerical  integrations 
of  these  data  and  compare  these  to  the  integrals  of  the  differential  equation 
representing  the  model  of  the  calorimeters.  For  the  backward  integrals  starting 
from  the  end  of  the  measurement  cycles  at  t  —  T,  we  obtain 

^  f^AH„et(r)dT  _  CpM[Ag(f)  -  A0(T)] 

""  SrM^)dr  Srfmdr 


while  the  corresponding  equation  for  forward  integration  from  the  start  of  the 
measurement  cycle  is 


_  AH„etiT)dr  Cr,M[Ae{t)  -  Ag(0)] 

lofiWdT  !lh{e)dT 


(14) 


Here,  the  suffices  21  and  31  denote  respectively  “backward  integration,  lower 
bound”  and  “forward  integration,  lower  bound.”  (fc^)2i  and  (fc^)3i  are  the  cor¬ 
responding  integral  heat  transfer  coefficients  defined  at  the  time  t. 


We  note  here  that  care  is  needed  when  integrating  the  term  [net  enthalpy  input, 
AjEf„et(r)]  around  the  discontinuities  a,t  t  =  ti  and  t  —  t2-  This  is  a  matter 
which  will  be  considered  further  in  Sections  5  and  6.  In  our  work  we  have  at 
various  times  used  the  trapezium  rule,  Simpson’s  rule,  or  the  mid-point  rule 
to  carry  out  the  integrations.  Of  these  rules,  only  the  mid-point  rule  is  strictly 

^While  dealing  with  these  “historical  aspects,”  we  note  also  that  in  1990-92  we  investigated 
the  use  of  Kalman  filtering  (leading  to  an  heat  transfer  coefficient  labelled  as  (A;^)6  and  we  also 
investigated  the  use  of  other  filtering  methods.  Some  of  this  was  rather  promising  especially 
that  designed  to  extract  information  about  “positive  feedback,”  but  these  projects  had  to  be 
abandoned.  As  part  of  these  projects  we  also  investigated  the  use  of  averaging  techniques 

other  than  the  “square-box”  version,  which  gives  from  and  from  (A;'^)ii. 

These  projects  were  also  promising,  but  again  had  to  be  abandoned.  Finally,  in  1994,  the 
change  of  the  ICARUS  systems  to  hardware  based  switching  was  investigated  with  a  view  to 
allowing  changes  in  the  data  acquisition  intervals  (thereby  putting  the  strategies  on 

a  sounder  basis).  However,  these  switching  systems  were  not  incorporated  in  the  ICARUS-2 
systems. 
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speaking  correct  in  that  it  agrees  with  the  mathematical  definition  of  an  integral. 
It  is  quite  generally  assumed  that  integrations  carried  out  using  the  trapezium  or 
Simpson’s  rule  will  converge  onto  the  “correct”  algebraic  result  if  the  integration 
interval  is  made  adequately  small,  but  this  does  not  necessarily  follow.  It  is  a 
matter  which  needs  to  be  investigated  for  each  particular  case. 


The  merits  of  the  particular  integration  procedures  coupled  to  the  adequacy  of 
the  chosen  integration  interval  is  revealed  more  clearly  when  we  come  to  the 
use  of  Eqs.  (13)  and  (14)  to  determine  CpM  and  to  carry  out  extrapolations  to 
remove  the  effects  of  the  second  term  on  the  R.H.S.  of  Eqs.  (13)  and  (14)  on  the 
corresponding  heat  transfer  coefficients.  We  will  follow  here  the  procedure  laid 
down  in  the  ICARUS-1  Handbook  [1]  where  the  integrations  were  restricted  to 
the  region  of  application  of  the  heater  calibration  pulse.  For  backward  integra¬ 
tion,  we  obtain 


flAH„et{T)dr  _  CpM[Ae{t)  -  Ae{t2)]  ,  „,,o^ 

ilumr  =  s:.hmr . ^ 

while  for  the  forward  integration  we  have 

AH„et{r)dT  CpM[Ae{t)  -  Ae{t2)  ,  „,,o^ 

— ft . .  7^ - = - rt  f 

j^Ji{e)dT  J^^fi{0)dT 


It  will  be  seen  in  Sections  5  and  6  (more  especially  in  Section  6)  that  Eq.  (15)  can 
be  used  to  derive  accurate  values  of  CpM  while  there  is  some  minor  degradation 
when  using  the  forward  integration,  Eq.  (16).  The  application  of  Eq.  (15) 
to  the  data  sets  was  the  “target  methodology”  of  the  ICARUS  systems  and 
the  derived  “lower  bound  heat  transfer  coefficient,”  (/:^^)26i  was  described  as 
in  the  Handbook  [1]  and  the  associated  correspondence.  The  same  types 
of  equation  may  be  used  to  derive  (A:^^)25i?  ik'^)27i,  and  281  as  well  as 
(^H^)37i7  ^nd  (A:if)38i;  it  is  only  necessary  to  start  the  interpretation 
from  the  appropriate  times,  which  also  give  the  starting  values  of  0  for  the  R.H.S. 
of  the  relevant  equations.  Of  these  sets  of  estimates,  that  leading  to  (/;:^^)28i  is 
especially  useful  and  this  particular  fit  also  gives  a  good  estimate  of  CpM. 


In  order  to  obtain  the  “true  heat  transfer  coefficients”  it  is  necessary  to  combine 
the  integrals  of  the  enthalpy  inputs  in  Eqs.  (15)  and  (16)  with  thermal  balances 
made  at  one  or  a  series  of  points.®  We  will  confine  attention  here  to  the  procedure 
originally  suggested  in  the  Handbook  for  the  ICARUS-1  system[l].  We  make  a 

®This  can  be  done  in  a  number  of  ways  and  it  is  important  that  this  part  of  the  evaluation 
was  changed  during  the  summer  of  1994  following  the  receipt  of  the  first  two  sets  of  data 
collected  by  NHE. 
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thermal  balance  just  before  the  application  of  the  calibration  pulse  and,  if  the 
system  has  relaxed  adequately  and,  if  dO/dt  =  0,  then  if  we  consider  {k'f^)32, 

0  =  [AH„M][t  -  h]  +  Qf[t  -  h]  -  ik'MiOb  +  -  et][t  -  h]  (17) 


Combination  with  Eq.  (14)  eliminates  the  unknown  rate  of  excess  enthalpy 
generation,  Q  / .  We  obtain 

/!/  \  •fti  ^^net{T')dT  —  [AHnet{tl)][i  “  ^l]  CpM[A9{t)  “  A0(ii)]  ,  . 

(fci?)32- - t  r .... - FT7^ - ‘  ^ 

j^j2{0)dr  Jt,h{0)dr 


The  corresponding  equation  for  (A;^)22  follows  from  Eq.  (18)  on  replacing  ti  by 

The  corresponding  equation  for  (A:^)22j  based  on  the  backward  integration  of 
the  data  sets,  follows  from  Eq.  (18)  on  replacing  ti  by  t2-  It  is  also  convenient 
to  rewrite  the  derived  equation  in  the  straight  line  form: 

^^net{'^)dT  —  [AHnet{t2)][i^  “  ^2]  _  CpM[A6{t)  —  A0{t2)]  ^ 

flf2{0)dr  Jlf2{e)dT  ^ 

(A:^^)262  Wets  the  version  of  (fc^)22j  which  we  used  in  our  investigations  prior  to 
the  construction  of  the  ICARUS-1  system.  As  we  did  not  want  to  discuss  the 
differences  between  these  two  versions,  we  also  labelled  {k*^)262  with  the  suffices 
22.  It  should  be  noted  that  Eq.  (19)  is  soundly  based  (in  a  mathematical  sense) 
in  that  the  extrapolation  to  [A6{t)  —  A6{t2)]  =  0  gives  the  value  of  (/c^^)262 
at  a  well  defined  time,  t  =  t2>  This  extrapolation  automatically  removes  the 
effect  of  the  term  CpM[{6{t)-d{t2)]l  f2{0)dr  on  the  heat  transfer  coefficient. 

This  was  one  objective  for  our  methodology  because  CpM  is  the,  least  accurate 
parameter  in  the  analysis;  the  application  of  Eq.  (19)  to  the  data  sets  was  the 
“target  methodology”  for  evaluating  the  “true  heat  transfer  coefficients.” 


While  it  is  also  possible  to  write  Eq.  (18)  in  the  form  (19)  to  give  (fc^^)362,  this 
method  of  analysis  is  not  useful  as  the  range  of  the  extrapolation  required  is  too 
long  as  will  be  shown  in  Sections  5  and  6,  For  this  reason  we  recommended  in 
the  Handbook  [1]  that  (A:^)32  be  evaluated  at  times  close  to  t  =  t2  using  Eq. 
(18).  However,  in  view  of  the  errors  in  the  determination  of  CpM ^  these  values 
of  (A:^)32  are  inevitably  less  accurate  than  those  of  {k^^)2e2- 

®We  note  that  NHE  did  not  follow  the  instruction  in  the  ICARUS- 1  Handbook  [1]  to  use 
2-day  measurement  cycles  and,  for  the  reduced  time  scales  of  1-day  cycles  in  particular,  it  is 
necessary  to  include  the  term  CpM{d/\$/dt)  in  the  thermal  balances,  Eq.  (17).  However,  NHE 
continued  to  use  the  original  form  of  the  equation.  They  also  did  not  follow  the  instruction  to 
evaluate  (fc'^)32  at  times  close  to  t2- 
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We  note  here  also  that  one  must  be  somewhat  careful  in  carrying  out  the  re¬ 
quired  linear  regression  fitting  procedures,  a  matter  which  is  considered  further 
in  Sections  5  and  6. 
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SECTION  3:  DIFFERENTIAL  EQUATIONS  GOVERNING  THE 
BEHAVIOR  OF  THE  CALORIMETERS:  SIMULATIONS  OF  THE 
TEMPERATURE-TIME  SERIES. 

It  has  been  established  that  at  low  to  intermediate  cell  temperatures  (say, 
30  <  0  <  80°C')  the  behavior  of  the  calorimeters  is  modelled  adequately  by 
the  differential  equation 

CpM{dM/dt)  =  [Ec{t)  -  Eth,t,]I  +  Qfit)  +  [AQH(t  -  h)  -  AQH{t  -  ta)] 

-(37/4F)  [{PIP*  -  P)]  [{Cp,D,o,g  -  Cp,D,o,i)Ae  +  L] 

-ik^nWt  [1  -  Jt]  [/i  {e)/et  +  4<AA0]  .  (20) 

With  the  calorimeters  used  in  the  ICARUS-type  investigations,  the  conductive 
contribution  to  heat  transfer  is  small.  We  have  therefore  assumed  that  this  term 
can  be  “lumped”  into  the  radiative  term  by  allowing  a  small  increase  in  the 
radiative  heat  transfer  coefficient: 

radiative  heat  transfer  =  (k'^)  [1  -  7<]  [(^6  +  A^)'*  —  0^]  .  (21) 

If  the  time  dependence  of  the  heat  transfer  coefficient  is  not  included  explicitly 
in  this  equation,  then 

radiative  heat  transfer  =  (k^)[(^b  +  A0)''  —  (22) 

where  the  radiative  heat  transfer  coefficient  (fe^)  now  shows  a  weak  time  depen¬ 
dence. 

In  calculating  the  rate  of  enthalpy  removal  by  the  gas  stream, 

(3//4F)[P/(P*  -  PmCp,D,o,g  -  Cp,D,o,i)Ae  +  L],  (23) 

we  have  always  assumed  that  the  partial  pressure  of  D2O  (or  H2O)  in  this 
gas  stream  can  be  calculated  using  the  Clausius-Clapeyron  equation  with  the 
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latent  heat  of  evaporation,  L,  being  that  at  the  boiling  point.  Evaporative 
cooling  only  becomes  a  major  term  at  temperatures  close  to  the  boiling  point 
(say,  at  A0  >  70^C)  where  these  two  assumptions  are  justified.  At  low  to 
intermediate  temperatures,  AHev{t)  is  a  minor  correction  term  so  that  errors 
due  to  the  two  assumptions  introduce  second  order  small  quantities  (for  a  further 
approximation,  see  below). 


In  carrying  out  simulations  to  be  analyzed  by  the  methods  outlined  in  Section 
2,  we  have  further  usually  assumed  that  the  rate  of  excess  enthalpy  generation 
is  zero, 

Qfit)  =  0,  (24) 

and  that  the  radiative  heat  transfer  term  can  be  linearized:  rate  of  heat  transfer, 

4{k'n)elM  =  iK)A0,  (25) 

i.e.,  that  the  heat  transfer  is  now  pseudo-conductive.  This  is  an  assumption 
which  we  also  used  (with  several  restrictions)  in  our  original  investigation  [8]. 
We  note  here  that  the  heat  transfer  coefficients  for  the  Dewar  cells  used  at  that 
time  were  up  to  twice  those  which  are  calculated  from  the  Stefan-Boltzmann 
coefficient  and  the  radiative  surface  area  so  that  we  had  to  assume  that  the 
conductive  contribution  was  appreciable.  We  attributed  this  conductive  contri¬ 
bution  to  conduction  across  the  nominal  vacuum  gap  due  to  inadequate  evacua¬ 
tion/baking  of  the  Dewars.  It  was  therefore  not  clear  whether  the  heat  transfer 
term  should  be  described  as  being  pseudo-radiative  or  pseudo-conductive  and 
the  experiments  had  to  be  carried  out  in  such  a  way  that  the  errors  due  to  the 
limiting  assumptions  were  below  those  of  the  experiment.  In  our  later  work 
(including  that  carried  out  with  the  ICARUS  systems),  we  ensured  that  the 
vacuum  in  the  Dewars  was  sufficiently  hard  so  that  the  radiative  heat  transfer 
coefficient  was  now  given  by  the  product  of  the  Stefan-Boltzmann  coefficient  and 
the  radiative  surface  area.  However,  it  is  evident  that  the  vacuum  in  some  of 
the  cells  used  in  the  NHE  investigations  had  become  rather  soft. 


The  correct  description  of  heat  transfer  from  the  cell  is  a  matter  which  requires 
further  investigation.  While  the  limiting  assumptions  introduce  small  errors,  we 
can  ensure  that  these  errors  are  less  than  those  due  to  the  experiment.  It  is 
therefore  better  to  use  these  assumptions  rather  than  to  attempt  to  separate 
this  term  into  the  radiative  and  conductive  contributions.  However,  a  better 
approach  might  well  be  that  we  should  calculate  the  radiative  term  and  then 
derive  the  conductive  contribution  from  the  calibrations  [3].  This  would  ensure 
that  we  do  not  introduce  an  additional  parameter  into  the  modeling  of  the  cells. 
It  should  be  noted  that  such  a  methodology  would  automatically  ensure  the 
linearization  of  Eq.  (20).  In  carrying  out  such  modified  procedures,  we  should 
weight  the  radiative  contribution  appropriately  (say  95%  of  the  total)  rather 
than  using  the  50:50  split  of  the  original  analysis  [3] . 
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The  use  of  the  descriptions  and  assumptions  outlined  above  gives  us  the  differ¬ 
ential  equation 

C^M^{dMldt)  =  [E,{t)  -  Eth,b]I  +  Qf{t)  +  ^QH{t  -  h)  --  AQH{t  -  ^2)- 

AH^t)  -  (4")[1  -  -  ^t]-  (26) 

In  order  to  calculate  the  Ad -t  series  we  have  to  deal  with  a  further  difficulty  in 
that  the  function  Ec{t)  is  unknown.  However,  it  is  evident  from  the  experimental 
A0(t)  -  t  and  Edt)  -  t  series  that  we  can  always  observe  “negative  feedback,” 
i.e.,  that  Ec  decreases  with  A6.  If  we  assume  that  this  is  the  only  cause  of  the 
variation  of  Ec  with  t,  then 

Ec{t)  -  Ec{0)  -  a{A9  -  -  E'  -  aAO  (27) 

(i.e.,  we  neglect  any  variation  of  the  activation  overpotentials  with  time)  and 
CpM{dA0/dt)  -  [£;;  -  Eth,b  -  otAO]!  +  Qf{t)  +  [AQH{t  -  h)  -  AQH{t  -  h)] 
-AHe.it)  -  (4°)(1  -  Jt)m  +  -  et].  (28) 


It  was  not  clear  during  1991-92  whether  we  should  change  our  data  processing 
strategy  to  that  later  incorporated  into  the  ICARUS  systems  or  whether  we 
should  continue  to  use  nonlinear  regression  fitting.  At  that  time,  we  therefore 
investigated  further  the  possibility  of  obtaining  an  analytic  solution  of  Eq.  (27) 
so  as  to  speed  up  the  latter  procedure.  Such  a  solution  was  derived  and  it  was 
shown  that  the  numerical  integration  of  Eq.  (28)  agreed  with  this  solution  to 
better  than  0.1%,  the  target  figure  for  our  data  analyses  (these  solutions  required 
the  assumption  Qf{t)  =  constant,  which  is  in  any  event  necessary  so  as  to  achieve 
calibrations  of  the  systems).  We  believe  that  the  residual  discrepancy  is  due  to 
deficiencies  in  the  numerical  techniques,  i.e.,  that  the  anal3d;ic  solutions  are  exact. 
This  program  of  work  was  discontinued  when  it  became  clear  that  satisfactory 
data  analysis  could  be  achieved  by  using  linear  regression  procedures.  However, 
it  may  well  be  that  this  particular  aspect  of  the  data  analysis  procedures  should 
be  restarted  if  we  wish  to  develop  a  general  investigation  of  existing  data  sets. 

At  that  time,  we  also  investigated  the  application  of  various  data  evaluation 
procedures  to  simulated  data  sets.  Attention  was  confined  to  linearized  versions 
of  Eq.  (28)  with  the  additional  assumptions, 


0 

11 

cy 

(29) 

AHevit)  =  0, 

(30) 

giving  the  equation. 
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CpMidM/dt)  =  E'J  +  [£^QH{t  -U)-  ^QH{t  -  fz)] 


(31) 


With  the  additional  assumption  that  the  pseudo-conductive  heat  transfer  coef¬ 
ficient  is  independent  of  time  we  have 

CpM{dA9/dt)  =  E'J  +  [AQH{t  -ti)-  AQHit  -  tz)]  -  (k'f  +  aI)A0.  (32) 


Parts  of  these  investigations  have  now  been  repeated.  The  integrations  of  (32) 
follow  immediately.  For  the  initial  condition 


o 

!l 

<1 

ll 

<1 

(33) 

we  obtain 

=  A6'ie--'‘ +  Ai(l  -  e-^‘) 

(34) 

where 

X  =  {k'J  +  aI)/CpM 

(35) 

and 

A^=E'J{k'f  +  aI) 

(36) 

i.e.,  in  the  region  0  <  t  <ti.  We  assume 

CpM°  =  330 J/f-i 

(37) 

k'f  =  0.073WK-^ 

(38) 

7/  =  QmiWK-^ 

(39) 

A9i  =  11.75“C' 

(40) 

E'l  =  JA9i  =  IW. 

(41) 

With  the  exception  of  the  rather  low  value  of  these  parameters  are  close 
to  those  which  we  would  derive  for  the  blank  experiment  discussed  in  Section  6. 
The  reason  why  A6i  has  been  set  artificially  low  is  so  as  to  allow  an  examination 
of  the  region  close  to  i  —  0.  With  (36)  through  (40),  we  obtain  from  Eq.  (34) 

A9  =  13.528125  -  1.778125et~°°°°(^'‘-^‘'  (42) 

while  the  enthalpy  input  is  given  by 

input  enthalpy  =  0.987553  -t-  0.012447exp[-0.000(24.)t].  (43) 
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Here  (24.)  denotes  the  recurrence  of  the  number  group  24. 


We  next  use  (41)  to  give  the  value  of  A^i  —  43,200s 

AOi  =  13.528075°C  at  t  =  43, 200s.  (44) 

This  value  is  adequately  close  to  13.528125®C,  which  applies  to  complete  thermal 
relaxation.  We  use  this  value  together  with 

AQ  =  0.2W,  ti  <t<t2  (45) 

for  the  next  step  of  the  integration.  With 

A2  =  {E'J  +  AQ)/ik'f  +  al)  (46) 

A9  =  +  A2[1  - 

=  16.028125  -  (47) 

and 

input  enthalpy  =  0.970053  +  0.017500e^”^‘^®^^^^'^^^”*^^^  (48) 

Equation  (47)  can  be  used  in  turn  to  derive 

Ai92  =  16.028049®C  at  (^2  -  ti)  =  42,900s,  (49) 

which  is  used  as  the  initial  value  for  the  final  step  in  the  integration.  We  obtain 

Ae  =  +  Ai[l  - 

=  13.528125  +  2.499924e>“°’°“(2^-)(‘"‘=^]  (50) 

and 

input  enthalpy  =  0.987553  -  0.017499el-°““(^'‘-^(*“*")].  (51) 


The  “raw  data”  calculated  using  Eqs.  (42),  (43),  (47),  (48),  (50),  and  (51)  are 
given  in  the  ICARUS  spreadsheet,  spreadsheet  1.  These  data  are  ana¬ 

lyzed  in  Section  5  using  the  methodologies  outlined  in  Section  3.  It  will  be  seen 
that  simulations  of  this  kind  are  adequate  for  demonstrating  the  advantages 
and  shortcomings  of  the  various  possible  methods.  It  will  also  be  clear  that 
such  simulations  are  deficient  in  several  important  respects.  In  the  first  place 
the  residual  time  dependence  of  the  heat  transfer  coefficient  has  not  been  taken 
into  account.  However,  it  was  confirmed  in  1992  that  a  data  set  generated  with 
approximate  solutions  which  include  this  effect  does  indeed  generate  the  key 
features  observed  when  using  “real  experimental  data,”  see  Section  6.  We  note 
here  that  data  sets  have  never  yet  been  generated  using  the  full  analytic  solution 
described  above.  Secondly,  the  effects  of  “noise”  having  defined  characteristics 
(i.e.,  defined  power  spectral  densities  as  well  as  the  effects  of  quantization  of  the 


21 


measurements)  on  the  methods  of  data  evaluation  has  never  yet  been  investi¬ 
gated.  Instead,  we  have  relied  on  a  comparison  of  the  analysis  of  simulated  data, 
Section  5,  with  those  for  “blank”  experiments.  Section  6.  If  the  question  of  the 
effects  of  “noise”  ever  becomes  an  important  issue,  then  we  would  suggest  that 
the  “noise”  characteristics  of  Ec{t)  and  A0(i)  be  first  determined  and  that  data 
be  then  generated  using  the  full  analytic  solution  with  addition  of  the  correct 
“noise”  characteristics. 
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SECTION  4:  SPECIFICATION  OF  THE  ICARUS-1  EXPERIMEN¬ 
TAL  PROTOCOLS  AND  DATA  EVALUATION  PROCEDURES. 


Before  dealing  with  the  analysis  of  the  simulated  data,  Section  5,  and  the  analysis 
of  data  for  a  “blank”  experiment,  Section  6,  we  will  outline  the  experimental  pro¬ 
tocols  and  data  evaluation  strategies  specified  for  experiments  with  the  ICARUS 
systems.  The  key  elements  were  as  follows: 


(i)  the  measurement  cycles  should  be  lengthened  to  48  hours.  Following  the 
replenishment  of  the  D2O  (or  H2O)  in  the  cells  to  make  up  for  losses  due  to 
electrolysis,  the  relaxation  of  the  systems  was  to  be  followed  for  12  hours  followed 
by  the  application  of  a  Joule  heating  pulse  for  a  further  12  hours  (when  this 
calibration  was  required)  in  turn  followed  by  a  final  relaxation  for  a  further 
24-hour  period. 

(ii)  the  protocols  were  to  be: 

(a)  two  measurement  cycles  without  calibration  pulses; 

(b)  ten  measurement  cycles  with  calibration  pulses; 

(c)  two  measurement  cycles  without  calibration  pulses; 

(d)  ten  measurement  cycles  with  calibration  pulses. 


It  will  be  seen  that  a  total  experiment  duration  of  48  days  was  specified; 


(iii)  several  “blank  experiments”  were  to  be  carried  out  (at  least  one  for  each 
cell  in  use;  the  use  of  Pt  cathodes  in  O.IM  Li0D/D2O  was  recommended).  The 
protocol  (ii)  was  to  be  followed, 

(iv)  the  execution  of  the  “blank  experiments”  was  to  be  followed  by  experiments 
using  cathodes  made  of  Johnson  Matthey  Material  Type  A.  The  protocol  (ii) 
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was  again  to  be  followed; 


(v)  the  first  step  in  the  data  evaluation  was  to  be  the  plotting  of  A3-  or  A4- 
sized  graphs  of  the  raw  data.  The  heat  transfer  coefficients  (A;^)i  and  (fc^)2 
were  to  be  derived  for  each  measurement  cycle. 


(vi)  the  next  step  was  to  be  the  construction  of  (fci^)ii-type  spreadsheets  cou¬ 
pled  to  the  determination  and  interpretation  of  and  (fc/^)ii.  The 

further  evaluation  of  these  spreadsheets  was  not  specified  in  1993;  this  was  a 
matter  which  was  to  be  decided  by  a  collaborative  program  between  NHE  and 
IMRA  Europe. 

(vii)  after  the  execution  of  (vi),  the  (A:^)2i-type  spreadsheets  were  to  be  pre¬ 
pared  and  values  of  (/u^^)26i  and  {k*^)36i  and  the  associated  values  of  CpM 
were  to  be  determined.  These  values  of  CpM  were  to  be  used  both  to  correct 
the  evaluations  in  (vi)  and  to  determine  the  “true”  heat  transfer  coefficients 
(A:^)32  at  times  close  to  the  end  of  the  calibration  period,  t  ~  t2\ 


(viii)  it  was  envisaged  that,  following  the  completion  of  this  initial  stage  of 
the  investigation,  the  ICARUS  program  would  move  on  to  the  examination  of 
materials  variables  as  well  as  the  production  of  ICARUS-2.  This  second  part  was 
intended  to  deal  with  the  effects  of  increasing  the  current  density,  the  analysis  of 
data  close  to  the  boiling  point  (including  the  boiling  episodes),  and  “Heat  after 
Death.” 


After  the  receipt  of  the  data  for  the  first  set  of  experiments  carried  out  in  the 
NHE  Laboratories,  it  became  apparent  that  there  were  timing  errors  in  the 
ICARUS-1  system  installed  in  Sapporo  [1],  The  most  self-evident  error  was  in 
the  timing  of  the  application  and  cessation  of  the  calibration  pulses  {ti  and  ^2), 
which  degraded  somewhat  the  estimation  of  (A:i^)3i  and  caused  a  serious  degra¬ 
dation  of  the  evaluation  of  (A:^)32-  The  estimations  of  ikj^)2i  and  {k'fi)22  were 
not  affected  by  these  errors  (as  had  been  the  case  for  the  experiments  carried 
out  prior  to  the  construction  of  the  ICARUS- 1  system).  It  was  therefore  recom¬ 
mended  [16]  that  (vii)  be  modified  and  that  the  “true  heat  transfer  coefficient” 
be  estimated  using  (A:^)22  evaluated  at  times  close  to  ti  and  by  (A;^^)262  (rather 
than  by  (A:^)32  estimated  at  times  close  to  t2  as  recommended  in  the  Handbook 
for  the  ICARUS-1  system  [1]), 

This  set  of  objectives  and  instructions  will  be  reconsidered  in  Section  7  in  the 
light  of  the  evaluation  of  the  simulated  data.  Section  5,  and  of  a  “blank”  exper¬ 
iment,  Section  6. 


24 


SECTION  5:  EVALUATION  OF  THE  “RAW  DATA”  GENERATED 
USING  THE  SIMULATION  DESCRIBED  IN  SECTION  4. 

(a)  We  start  by  plotting  the  “raw  data”  generated  by  the  simulation  Eqs. 
(42),  (43),  (47),  (48),  (50),  and  (51)  described  in  Section  3;  see  the  ICARUS 
(fc')ii “Spreadsheet,  spreadsheet  1.  The  plots  of  the  temperature-time  and  input 
enthalpy-time  series  are  shown  in  Fig.  2  (the  data  are  the  same  as  in  Fig.  1). 
We  evaluate  (fc')i  and  (A;')2  using  the  graphical  method  and  obtain 

(fc;)i  -  0.07294WA:"^  (52) 

(fc;)2  =  0.07280 (53) 


It  should  be  noted  that  the  scale  of  the  y-axis  in  Fig.  2  is  markedly  reduced 
compared  to  that  used  in  the  evaluation  of  experimental  data,  e.g.,  Fig.  38. 
The  reasons  for  this  are  the  rather  low  value  of  A6i  used  in  the  simulation  (Eq. 
(39))  which  increases  the  excursion  of  the  temperature-time  series  coupled  to 
use  of  landscape  rather  than  portrait  format  so  as  to  allow  the  presentation 
of  the  whole  of  the  measurement  cycle.  In  view  of  this  compressed  scale,  the 
accuracies  of  and  (fcc)2  are  somewhat  reduced  compared  to  those  normally 
achieved.  This  effect  is  to  some  extent  counteracted  by  the  absence  of  “noise” 
in  the  simulated  data. 

At  this  stage,  we  have  also  always  prepared  a  diagram  on  a  much  larger  scale  of 
the  regions  straddling  the  times  ti  and  t2  such  as  that  shown  in  Fig.  3.  There 
are  several  reasons  for  preparing  such  diagrams.  First  of  all,  when  dealing  with 
experimental  data  collected  with  the  ICARUS  systems  (and  their  precursors), 
it  is  important  to  determine  the  exact  times  of  application  and  cessation  of  the 
calibration  pulses.  It  was  the  preparation  of  such  diagrams  which  showed  that 
there  were  timing  errors  in  the  ICARUS- 1  system  which  could  be  easily  allowed 
for  by  deriving  the  exact  values  of  ti  and  ^2-  As  has  already  been  stated  on 
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other  occasions  there  was  no  point  in  changing  the  hardware  to  correct  for  these 
effects  as  the  times  of  data  acquisition  were  known  with  sufficient  accuracy. 

Of  course,  when  dealing  with  the  simulated  data,  ti  and  t2  are  known  exactly  and 
these  times  should  also  be  known  accurately  for  experimental  data  collected  with 
the  ICARUS-2  system.  However,  we  believe  that  this  assumption/prediction  will 
turn  out  to  be  incorrect! 


The  second  reason  for  investigating  the  temperature-time  series  in  these  time 
regions  is  to  define  the  time  series  required  for  the  integrations  of  the  input 
enthalpy  to  be  used  in  the  (fc^)2i-spreadsheets.  This  definition  is  in  turn  de¬ 
pendent  on  the  methods  to  be  used  for  the  numerical  integrations.  Here,  we  will 
assume  here  that  we  have  decided  to  use  the  trapezium  or  sum  rules  for  the  in¬ 
tegrations  (these  methods  are  related  because  the  values  given  by  the  trapezium 
rule  are  those  derived  by  the  sum  rule  less  one  half  the  terms  at  the  extremes 
of  the  range  multiplied  by  the  integration  interval).  For  these  rules,  we  need  to 
insert  an  additional  point  into  the  experimental  data  at  and  t2  if  these  times 
correspond  to  the  data  acquisition  intervals.  If  ti  and  t2  do  not  correspond  to 
these  intervals,  we  need  to  insert  two  additional  points  at  each  of  ti  and  t2  • 


We  note  here  that  the  trapezium  and  sum  rules  can  be  used  in  a  straightforward 
way  to  integrate  around  the  discontinuities  at  ti  and  ^2-  The  mid-point  and 
Simpson’s  rule  cannot  be  used  in  this  way:  in  particular,  it  is  necessary  to 
use  much  more  complicated  procedures  if  we  wish  to  use  the  mathematically 
exact  central  difference  methods.  For  this  reason,  we  have  always  relied  on 
the  trapezium  or  sum  rules  except  if  the  integrations  are  confined  to  regions 
0  <  ^  <  i  <  ^2  or  ^2  <  ^  < 

(b)  The  next  step  is  to  carry  out  a  detailed  examination  of  the  {k'^)n-  spread¬ 
sheet,  spreadsheet  1  and  this  examination  itself  falls  into  several  parts.  We  first 
of  all  plot  the  lower  bound  heat  transfer  coefficient,  (A:c)ii5  against  time  (here 
given  by  the  measurement  interval)  where  we  use  the  assumed  value  of  CpM, 
Fig.  4.  Of  course,  for  the  case  of  the  interpretation  of  the  simulated  data,  we 
know  the  value  of  CpM  used  in  the  simulation.  The  deviation  of  (A;^)ii  from 
the  value  0.073000  WK“^  used  in  the  simulation  is  therefore  an  indication  of 

believe  that  these  integrations  have  been  carried  out  incorrectly  by  the  ICARUS-2 
software.  It  would  be  important  to  check  this  particular  point  if  we  can  reimplement  the 
software.  The  reason  is  not  only  that  incorrect  integrations  will  lead  to  incorrect  values  of 
the  “true  heat  transfer  coefficient,”  a  matter  which  will  be  illustrated  below,  but  also  that  we 
will  calculate  an  incorrect  value  of  the  excess  enthalpy  for  each  measurement  cycle  even  if  we 
should  have  a  correct  value  for  the  heat  transfer  coefficient!  If  we  cannot  reimplement  the 
software,  then  we  will  be  able  to  make  some  sort  of  overall  check  of  the  ICARUS-2  software 
by  carrying  out  correct  and  detailed  checks  of  the  measurement  cycles  for  which  NHE  have 
given  values  of  the  “true  heat  transfer  coefficient”  evaluated  by  their  methodology.  This  should 
automatically  reveal  the  method(s)  they  actually  used  to  arrive  at  their  conclusions. 
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the  errors  in  dAd/dt  produced  by  using  the  second  order  central  difference  when 
using  a  300  s  measurement  interval.  We  can  see  that  we  can  obtain  reasonable 
values  of  (jfcc)ii  provided  we  exclude,  say,  the  first  six  hours  following  any  change 
in  the  operating  conditions  of  the  cells. 


Of  course,  in  the  case  of  the  evaluation  of  spreadsheets  derived  for  experimental 
measurements,  our  first  evaluation  will  be  based  on  a  “guesstimate”  of  CpM. 
It  will  therefore  be  necessary  to  amend  the  spreadsheet  once  an  accurate  value 
of  CpM  has  been  derived.^^  It  is  important  that  under  such  conditions  the 

further  average  (fc^)ii  cannot  be  derived  and  (A:')ii  is  quite  meaningless.  It 
is  therefore  impossible  to  carry  out  these  important  averaging  methods  which 
reduce  the  effects  of  random  errors.  These  difficulties  could  have  been  avoided  if 
the  experiments  had  been  carried  out  using  the  stipulated  48-hour  measurement 
cycles. 

We  note  that  the  error  of  0.04%  in  the  “lower  bound  heat  transfer  coefficient”  is 
above  that  which  was  specified  for  the  ICARUS-1  system  (0.01%).  An  alternative 
approach  to  removing  the  errors  due  to  incorrect  estimates  of  CpM  is  to  use  the 
evaluations  of  {k'f)i5u  (fcc°)i6i,  k'f)i7u  and  Figs.  5  through  8.  These 

extrapolations  automatically  give  us  values  of  CpM  and  the  relevant  values  of 
the  “lower  bound  heat  transfer  coefficients”  are  shown  on  Fig.  4. 

This  particular  approach  was  de-emphasized  in  setting  up  the  IC  ARUS-l  system 
because  the  extrapolations  in  Figs.  5  trough  8  are  to  a  point  where  dA0/dt  —  0. 
While  this  condition  is  satisfied  for  the  simulations,  it  will  not  be  observed  for 
experiments  carried  out  with  a  24-hour  measurement  cycle.  If  a  48-hour  cycle  is 
used,  then  there  will  be  two  times  at  which  this  condition  will  hold  (one  within 
the  period  0  <  ^  <  and  one  within  h  <  t  <  see  Section  6).  However,  it 
has  never  been  established  that  the  plots  in  Figs.  5,  6  and  8  extrapolate  to  these 
points.  This  is  a  matter  which  should  be  investigated  using  data  generated  by 
more  elaborate  simulations. 


(c)  The  next  step  is  to  investigate  the  “true  heat  transfer  coefficient”,  (fc^)i2. 
In  the  region  ti  <  i  <  ^2,  we  can  apply  Eq.  (9)  at  any  chosen  point.  This  is 
straightforward  for  the  data  derived  by  simulation  because  the  thermal  balance 
in  the  absence  of  the  calibration  pulse  is  known  exactly.  For  experimental  data, 
the  relevant  values  can  be  obtained  from  the  plots  of  the  temperature-time  and 
cell  potential-time  series  such  as  those  shown  in  Fig.  38,  Section  6.  Figure  9 
gives  the  relevant  values  using  the  data  on  the  (/c')ii-spreadsheet,  spreadsheet  1. 

^^The  little  which  we  have  seen  of  the  NHE  evaluations  leads  us  to  think  that  this  further 
step  was  never  carried  out.  If  this  is  correct,  then  the  values  of  CpM  used  would  have  been  in 
error  by  between  10  and  20%.  This  error  would  in  turn  have  produced  an  error  of  -  0.04%  in 
the  “lower  bound  heat  transfer  coefficients”  under  the  most  favorable  conditions  of  estimating 
this  coefficient  at  the  end  of  the  relevant  6-hour  periods. 
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The  deviations  from  the  known  value  ((A:c)i2  ==  0.073000  WK“^)  shown  in  Fig.  9 
as  well  as  the  component  parts  of  this  heat  transfer  coefficient  (the  ordinates  and 
abscissae  shown  in  spreadsheet  1)  show  that  the  effects  of  errors  in  the  estimates 
of  {d^d/dt)  are  now  more  serious  than  is  the  case  for  the  calculation  of  (fcc)ii* 
This  is  to  be  expected  mainly  in  view  of  the  relative  magnitudes  of  /i  [6)  and 
/2(0),  Eqs.  (9)  and  (10).  However,  we  can  see  from  spreadsheet  2  that  we  would 
be  able  to  derive  values  of  (fc' )i2  accurate  to  0.01%  if  measurement  cycles  lading 
48  hours  were  used.  It  is  possible  to  apply  averaging  procedures  leading  to  (fcc)i2 
and  (^)i2  for  the  last  6  hours  of  calibration  periods  lasting  12  hours,  thereby, 
markedly  reducing  the  effects  of  “noise”  in  the  experimental  temperature-  and 
cell  potential-time  series.  The  benefits  of  this  type  of  averaging  are,  indeed, 
foreshadowed  by  the  accurate  determination  of  Fig.  2. 


The  restrictions  which  should  be  placed  on  the  determination  of  the  “differen¬ 
tial  true  heat  transfer  coefficients”  are  also  illustrated  by  the  determination  of 
(A:'’^)i62,  Fig.  10.  It  can  be  seen  from  this  figure  and/or  spreadsheet  2  that  the 
range  of  the  extrapolation  required  to  remove  (nominally)  the  effects  CpM 
on  the  derived  heat  transfer  coefficient  is  much  longer  than  for  the  equivalent 
determination  of  (A:^’^)i6i.  The  abscissae  are  still  ca  10%  of  the  ordinates  even 
at  the  end  of  the  region  which  can  sensibly  be  used  for  the  determination  of 
(A:'’")i62  (say,  3  hours  after  the  application  of  the  calibration  pulse).  It  follows 
that  a  20%  error  in  CpM  will  lead  to  at  least  a  2%  error  in  (fc'’®)i62*  Of  course, 
such  effects  are  not  apparent  when  using  data  produced  by  simulations  free  of 
“noise”  but  we  would  predict  that  procedures  based  on  making  thermal  balances 
close  to  t  will  not  give  accurate  “true  heat  transfer  coefficients”  and  that 
we  should  instead  rely  on  evaluating  {k[)i2  (and  of  related  coefficients)  at  times 
close  to  i  —  ^2-  This  point,  which  will  be  further  illustrated  in  Section  6,  was  of 
key  importance  to  the  specification  of  the  ICARUS  procedures  [1]. 

The  value  of  obtained  from  Fig.  10  is  also  shown  on  Fig.  9.  Although 

this  value  agrees  with  those  described  by  methods  judged  to  be  satisfactory,  it 
should  be  noted  that  this  agreement  is  largely  due  to  the  use  of  data  free  from  all 
ambiguities  (absence  of  “noise,”  time  dependence  of  the  heat  transfer  coefficients, 
timing  errors).  The  use  of  experimental  data  leads  to  a  marked  degradation  of 
the  evaluation,  see  Section  6.  A  major  reason  for  this  degradation  is  that  the 
data  points  having  the  highest  statistical  weight  are  also  the  ones  which  have 
the  lowest  accuracy.  Fig.  10  and  spreadsheet  1.  It  is  possible,  however,  that 
the  use  of  statistical  weighting  procedures  would  allow  satisfactory  estimates  of 
(A:'’^)i62  to  be  made. 

Figure  9  also  shows  the  values  of  (A:'’^)i52  and  (fcc’^)i72  derived  from  spreadsheet 
1  by  using  Eq.  (12).  The  form  of  this  equation  shows  that  the  derivation  of  these 
versions  of  the  heat  transfer  coefficient  is,  in  fact,  dependent  on  the  interpretation 
of  cooling  (or  heating)  curves,  i.e.,  the  driving  force  is  now  the  enthalpy  content 
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of  the  calorimeters  rather  than  the  heater  calibration  pulse,  AQ  (these  methods 
were  developed  largely  as  part  of  a  program  of  work  on  the  interpretation  of 
phenomena  linked  to  “Heat  after  Death”).  As  would  be  expected,  the  accuracy 
of  these  estimates  of  the  heat  transfer  coefficient  is  limited  even  when  using  data 
free  from  all  ambiguities.  The  accuracy  of  estimates  of  {kc’^)i52  and  (A:c’^)i62 
can  never  be  greater  than  that  of  CpM  when  using  real  experimental  data.  This 
accuracy  is  further  degraded  by  the  fact  that  the  estimates  are  dependent  on  the 
derivation  of  {dA6/dt)  in  regions  where  the  accuracy  of  this  gradient  is  limited 
by  the  length  of  the  data  acquisition  interval. 

Although  Fig.  9  shows  three  data  points  for  each  of  these  estimates  of  the 
“true  heat  transfer  coefficient,”  we  would  expect  that  only  one  such  data  point 
could  be  derived  when  using  experimental  data  in  view  of  the  need  to  estimate 
{dA6/dt). 


d)  In  this  sub-section,  we  are  also  including  a  second  set  of  evaluations  which 
exactly  follow  the  set  described  above,  Figs.  4  through  10,  except  that  dAO/dt 
has  been  estimated  using  the  first  order  backward  difference.  Comparison  of 
Figs.  11  and  4  shows  that  the  errors  in  are  now  markedly  increased  in 

regions  where  {dAO/dt)  is  large:  as  expected,  the  errors  in  (dAO/dt)  based  on 
the  first  order  differences  are  much  larger  than  those  based  on  the  second  order 
differences.  Nevertheless,  satisfactory  values  of  can  be  obtained  provided 

measurement  cycles  of  48-hour  durations  are  used. 

Figure  11  also  includes  values  of  fcc^)i6ij  fcc^)i7i  and  (/c'’°)i8i  deter¬ 

mined  by  using  the  relevant  extrapolations,  Fig.  12  through  15.  It  can  be  seen 
that  these  values  of  the  heat  transfer  coefficient  can  again  be  estimated  satis¬ 
factorily.  However,  the  slopes  of  the  regression  lines  are  markedly  reduced  (due 
to  the  incorrect  estimation  of  (dAO/dt))  so  that  the  values  of  CpM  are  now 
also  too  low.  Of  course,  application  of  this  methodology  to  “real  experimental 
data”  would  then  lead  to  erroneous  estimates  of  (/c')ii-  Figure  16  gives  the  plot 
of  (fc')  12  versus  time  again  based  on  the  first  order  backward  differences  while 
Fig.  17  illustrates  the  estimation  of  (A;'’^)i62-  Comparison  of  Figs.  16  and  9 
shows  that  (k'^)i2  is  now  markedly  in  error  in  the  time  region  adjacent  to 
it  is  necessary  to  extend  measurements  to  at  least  7  hours  in  order  to  ensure 
that  the  errors  due  to  the  incorrect  estimation  of  (dAO/dt)  fall  below  0.1%  (the 
target  specification  for  ICARUS-1).  Of  course,  this  situation  if  aggravated  when 
using  “real  experimental  data”  in  view  of  errors  in  estimation  of  CpM.  Figure  17 
shows  that  satisfactory  values  of  (fcc’^)i62  can  still  be  obtained  when  using  such 
incorrect  estimates  of  (dAO/dt).  This  would  be  expected  because  the  extrapo¬ 
lation  procedure  removes  the  effects  oi  CpM.  Needless  to  say,  the  value  of  CpM 
derived  is  markedly  in  error.  Comparison  of  Fig.  17  with  spreadsheet  2  leads  to 
a  further  important  conclusion.  We  can  see  that  satisfactory  extrapolations  can 
be  obtained  even  though  the  “point-by-point”  values  are  totally  in  error  at  short 
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times.  As  has  already  been  noted,  such  long  extrapolations  are  to  be  avoided 
when  evaluating  experimental  data. 

Figure  16  also  includes  some  values  of  (A:c’^)i52  and  (A:'’^)i72  derived  from  spread¬ 
sheet  2.  It  can  be  seen  that  these  values  are  markedly  in  error  even  when  using 
the  correct  value  of  CpM.  However,  we  note  that  values  close  to  0.073000  WK“^ 
(the  value  in  the  simulation)  would  be  obtained  if  the  incorrect  estimate  of  the 
water  equivalent  (318.1  JK“^)  derived  in  Fig.  17  were  used  to  calculate  the 
“point-by-point”  values  of  (fcc’^)i52  and  (A:'’^c)i72.  It  follows  that  the  estimates 
of  CpM  and  of  (fcc’^)i52,  (fcc^)i62  and  (fcc’^)i72  are  internally  consistent  even 
through  the  component  parts  deviate  from  the  true  values.  The  procedures  used 
to  calculate  these  values  of  the  “true  heat  transfer  coefficients”  must  therefore 
be  used  with  due  care  -  in  fact  it  is  best  to  avoid  such  methodologies. 

(e)  We  will  consider  next  the  derivation  of  the  various  versions  of  the  integral 
heat  transfer  coefficients.  The  (A:^’°)2i  spreadsheet,  spreadsheet  3,  gives  the 
integrals  required  for  the  evaluation  of  and  (/!:')3i  and  Fig.  13  gives  the 

plots  of  these  two  heat  transfer  coefficients  against  the  time.  It  can  be  seen  that 
except  for  small  deviations  of  (A;')3i  from  the  value  0.073000  WK”^  used  in  the 
calculation  of  the  “raw  data,”  the  two  estimates  of  the  heat  transfer  coefficient 
agree  with  this  value.  The  reason  for  the  small  deviations  of  (A:')2i  at  short 
times  are  immediately  evident.  The  term 

CpM{e-0o)/  ['  f{e)dr 
Jt2 

is  negligibly  small  for  the  estimation  of  (A:^)2i  because  {6  —  ^o)  is  itself  small  for 
the  backward  integration.  On  the  other  hand,  the  term 

CpM{e  -  9o)/  f  fi0)dT 
Jtl 

is  initially  more  than  10%  of  (A:' )3i  for  the  forward  integration  procedure.  The 
deficiencies  of  using  the  trapezium  rule  coupled  to  the  use  of  an  inadequately 
long  data  acquisition  interval  are  therefore  immediately  apparent.  It  should  be 
noted  that  the  plots  of  (A:')2i  and  (A:')3i  versus  time  do  not  show  any  effect 
due  to  the  discontinuities  at  t  =  ti  and  t  =  provided  the  integral  of  the 
enthalpy  input  has  been  correctly  estimated  at  these  points.  The  reason  for  this 
suppression  of  the  effects  of  the  discontinuities  is  simply  that  the  magnitudes 
of  the  integrals  are  now  sufficiently  large  so  that  errors  due  to  the  use  of  the 
trapezium  rule  coupled  to  the  use  of  an  inadequate  long  measurement  interval 
are  no  longer  detectable. 

These  effects  (and  noneffects!)  due  to  errors  in  the  estimation  of  the  integrals 
were  of  key  importance  to  the  evolution  of  the  ICARUS- 1  data  processing  strat- 
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for  (*'’0)3 J, 2- 


Figure  19  and  spreadsheet  4  give  comparisons  of  (fcc)2i  and  (fcc)22-  It  can  be  seen 
that  the  evaluation  of  the  “true  heat  transfer  coefficient,  {k'j22,  is  satisfactory 
under  all  conditions.  This  is  also  brought  out  very  clearly  by  considering  the 
extrapolation  procedures,  Figs.  20  and  21.  In  particular,  in  the  estimation  of 
(*c’°)262,  the  term  ^ 

CpM{$-9o)/  [  f2iO)dT 
Jt2 

never  exceeds  10%  of  the  term 

f  AHdr/  [  f2{e)dT. 

Jt2 


The  situation  is  radically  different  for  the  evaluation  of  (fc'’^)362,  spreadsheet 
5.  While  the  extrapolation  procedure  for  obtaining  Fig  22,  is  still 

reasonably  satisfactory,  provided  the  points  immediately  adjacent  to  ti  are  ex¬ 
cluded  (the  abscissae  never  exceed  15%  of  the  ordinates,  spreadsheet  5),  that 
for  (fc^’°)362>  Fig.  23,  is  unsatisfactory  because  the  range  of  the  extrapolation 
required  is  very  long  (note  the  magnitudes  of  the  ordinates  and  abscissae  in 
spreadsheet  5  and  the  fact  that  the  abscissae  are  almost  equal  to  the  ordinates 
at  short  times).  Such  evaluations  will  also  be  markedly  degraded  by  the  “noise” 
of  real  experimental  data.  Furthermore,  even  the  borderline  fits  obtained  must 
be  to  some  extent  fortuitous  because  the  replacement  of  integrations  using  the 
trapezium  rule  by  the  mathematically  sound  mid-point  rule.  Figs.  24  and  es¬ 
pecially  Fig.  25,  gives  less  satisfactory  results  than  those  obtained  in  Fig.  23. 
(Note  especially  the  erroneous  values  of  CpM.) 


The  conflicting  effects  of  changes  in  the  integration  methods  and  data  acquisition 
interval  have  never  yet  been  resolved.  In  view  of  this  situation,  it  was  recom¬ 
mended  in  the  Handbook  for  the  ICARUS-1  system  [1]  that  (fcc)32  be  evaluated 
at  times  close  to  t2  where  the  effects  of  the  term 


CpM{e-eo)/  f  f2{e)dr 

Ju 


on  the  overall  value  of(fcJ.)32  are  reduced.  However,  following  the  receipt  of 
the  first  set  of  data  from  the  NHE  Laboratory  [11],  it  became  clear  that  the 
Group  would  never  achieve  satisfactory  evaluations  of  this  version  of  the  “true 
heat  transfer  coefficient.”  As  (fc'’^)262  could  be  evaluated  satisfactorily  under  all 
conditions,  it  was  recommended  that  future  evaluations  should  be  based  on  this 
version  of  the  “true  heat  transfer  coefficient.” 
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This  recommendation  was  also  influenced  by  a  consideration  of  the  effects  of 
timing  errors.  As  has  already  been  explained,  the  times  of  application  and 
cessation  of  the  heater  calibration  pulses  did  not  coincide  exactly  with  the  data 
acquisition  points  in  the  pre-ICARUS  phase  of  the  investigation.  Furthermore, 
there  were  evident  timing  errors  in  the  ICARUS-1  systems.  These  timing  errors 
did  not  affect  the  determination  of  (fc^’®)262  as  the  times  of  the  data  acquisition 
points  were  known  exactly.  Nevertheless  we  investigated  the  effects  of  gross 
errors  in  the  estimation  of  f  AH dr  on  the  evaluation  of  (fc'’^)26i  and  (fc'’^)262- 
See  Figs.  26  through  29  and  spreadsheet  4,  It  can  be  seen  that  such  gross  errors 
lead  to  a  maximum  error  of  0.5%  in  (/?c’^)262- 

Such  evaluations  are  therefore  reasonably  satisfactory.  On  the  other  hand,  eval¬ 
uations  of  (A:^’°)36i  and  (fc'’^c)362  fail  almost  completely.  Thus,  Figs.  30  and  31 
(spreadsheets  5)  show  that  erroneous  values  of  both  the  “lower  bound  heat  trans¬ 
fer  coefficients”  and  of  the  water  equivalent  are  obtained  if  there  are  errors  in 
f  AHdr.  Nevertheless,  the  values  of  (A:c)3i  at  long  times  are  still  within  the  ac¬ 
ceptable  range  as  specified  in  the  ICARUS-1  Handbook  [1].  The  extrapolations 
required  for  (A:'’^)362  fail  completely  (and  the  estimates  of  CpM  vary  widely)  as 
shown  by  Figs.  32  through  37,  Examination  of  the  sets  Figs.  32  through  34  and 
Figs.  35  through  37  shows  that  the  values  of  (A:'’^)362  and  CpM  deduced  depend 
on  the  range  of  the  regression  lines  fitted.  This  is  just  the  sort  of  behavior  which 
was  detected  in  the  limited  information  available  about  the  experiments  carried 
out  by  NHE,  We  believe  therefore  that  the  malfunctions  which  they  have  re¬ 
ported  are  in  large  measure  due  to  the  combination  of  estimating  the  “true  heat 
transfer  coefficients”  by  using  (fc^^)362  without  any  correction  of  the  effects  due 
to  timing  errors. 
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SECTION  6:  EVALUATION  OF  A  MEASUREMENT  CYCLE  FOR 
A  “BLANK  EXPERIMENT”  (Pt  cathode  polarized  in  O.IM  LiOD/DsO) 
USING  AN  ICARUS-2  SYSTEM. 

It  has  been  explained  elsewhere  the  reasons  why  we  have  very  few  data  sets 
collected  with  the  ICARUS-1  and  -2  systems.  However,  some  of  the  “raw  data” 
for  parts  of  a  “blank  experiment”  carried  out  during  the  summer  of  1995  using 
an  ICARUS-2  system  installed  at  IMRA  Europe  are  available.  This  experiment 
used  a  Pt  cathode  polarized  in  0.1  M  Li0D/D20  in  an  ICARUS-2  cell  (having  an 
extended  length  of  silvering  in  the  upper  portion  of  the  cell).  This  measurement 
cycle  belongs  to  one  set  of  calibrations  of  the  ICARUS-2  systems  carried  out  at 
that  time.  Unfortunately,  these  calibrations  were  terminated  in  1995  and  most 
of  the  data  collected  at  that  time  are  no  longer  available. 

This  data  set  can  be  regarded  as  being  satisfactory  except  in  one  regard.  At  that 
time,  the  cells  had  been  wired  to  the  “switching  boxes”  using  thin  wire.  This 
was  to  be  replaced  by  thick  wire,  but,  unfortunately,  the  wiring  of  the  calibration 
heaters  was  not  changed.  The  power  delivered  to  these  heaters  therefore  has  to 
be  corrected  for  voltage  losses  external  to  the  cell.  This  correction  is  ca  2%  of 
the  nominal  power  delivered  to  the  cell  and  can  be  made  to  better  than  0.1%  of 
the  2%  level.  The  possible  error  in  the  calibration  power  is  therefore  well  below 
the  target  value  for  obtaining  the  “true  heat  transfer  coefficients”  with  errors 
below  0.1%.  The  power  delivered  to  the  calibration  heater  was  0. 23025 W  and 
the  cell  current  was  0.20306 A. 


(a)  We  again  start  by  plotting  the  “raw  data,”  this  time  of  the  cell  temperature 
and  cell  potential  versus  time,  Fig.  38,  and  by  constructing  the  relevant 
spreadsheet,  spreadsheet  6.  Here,  we  have  used  the  third  measurement  cycle  so 
as  to  be  consistent  with  later  evaluations  (to  follow).  It  can  be  seen  that  such 
measurements  for  “blank  systems”  do  not  show  any  anomalies.  In  particular. 
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the  initial  temperature  perturbation  due  to  the  refilling  of  the  cell  (to  make 
up  for  the  losses  of  D2O  due  to  electrolysis)  relaxes  within,  say,  7  hours;  in 
view  of  the  much  larger  amplitude  of  the  temperature  perturbation  due  to  the 
application  and  cessation  of  the  heater  calibration  pulse,  we  need  to  allow  at 
least  8.5  hours  to  “eye-ball”  the  relaxations  in  the  regions  ti  <  t  <  t2  and 
t2  <  t  <T,  One  reason  for  the  specification  of  the  48-hour  long  measurement 
cycles  will  be  apparent  immediately;  such  measurement  cycles  allow  us  to  make 
the  durations  of  ti  ~  to  and  ^2  “  ^1  equal  to  12  hours. 


It  will  also  be  evident  that  the  temperature  following  the  cessation  of  the  cali¬ 
bration  pulse  relaxes  to  the  sloping  base  line.  Furthermore,  it  will  be  clear  that 
the  system  normally  shows  “negative  feedback”  in  that  increases  of  the  cell  tem¬ 
perature  lead  to  a  lowering  of  the  enthalpy  input  (and  vice  versa).  Tests  of  this 
“negative  feedback”  will  be  discussed  in  later  sections. 

If  we  accept  data  sets  such  as  those  in  Fig.  38  as  being  “reasonably  normal,”  we 
can  evaluate  the  “lower  bound”  and  “true  heat  transfer  coefficients”  at  a  time 
close  to  t2  by  using  the  graphical  methods.  We  obtain 

(A:^)ii  -  0,61844  x  with  Eth,b  -  1-54^  (54) 

(jfc^^)!!  =  0.62006  X  IQ-^WK-'^  with  Eth,b  =  1.527V  (55) 

(fc^)2  =  0.62027  X  10"^WK”^  (56) 


The  close  agreement  of  Eqs.  (54)  and  (56)  is  almost  certainly  fortuitous  because 
we  expect  the  errors  in  {k'ji)2  to  be  ca  ±  0,002  x  10“^  WK"^.  If  we  accept  the 
values  given  by  Eqs.  (55)  and  (56),  then  we  deduce  a  rate  of  excess  enthalpy 
generation  of  0.00034W,  a  value  which  is  comparable  to  those  which  we  had 
observed  previously  for  “blank  experiments”  [5,  6,  8].  As  was  noted  previously, 
such  low  values  are  below  those  expected  for  the  reduction  of  electrogenerated 
oxygen  present  in  the  solution.  This  can  be  explained  by  the  degassing  of  the 
solution  adjacent  to  the  cathode  by  the  electrogenerated  deuterium. 


In  this  preliminary  assessment,  we  also  prepare  plots  on  an  expanded  scale  of  the 
temperature-time  data  in  the  regions  adjacent  to  ti  and  ^2?  Fig.  39  (compare 
Fig.  3).  We  can  see  that  the  times  of  application  and  cessation  of  the  calibration 

As  we  have  always  explained,  it  is  essential  to  produce  such  graphs  of  the  raw  data  in  order 
to  check  on  the  normality  (or  otherwise)  of  the  experiments.  The  data  shown  in  Fig.  38  are 
“reasonably  normal,”  i.e.,  we  can  judge  them  as  being  suitable  for  further  evaluation  although 
it  was  apparent  that  the  “noise  levels”  had  increased  compared  to  the  data  sets  which  had 
been  collected  in  Salt  Lake  City  and,  subsequently,  in  the  “old  part”  of  the  IMRA,  Europe 
Building  in  Sophia  Antipolis.  This  increase  in  the  “noise  levels”  caused  some  degradation  in 
the  evaluations  compared  to  those  which  we  could  achieve  previously  up  to  the  end  of  1992. 
This  was  already  apparent  in  1994  at  which  time  we  tried  to  find  the  causes  for  the  increases 
in  “noise.”  Unfortunately,  the  attempts  to  do  so  were  stopped. 
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pulse  are  adequately  synchronized  with  the  data  acquisition  points  (although 
such  extrapolations  for  “real”  experimental  data  are  inevitably  less  certain  than 
those  for  simulation,  Fig.  3).  It  follows  that  it  is  sufficient  to  insert  single 
additional  points  at  h  and  t2  in  the  preparation  of  the  (A;^)2i -spreadsheet. 


The  evaluation  of  (A:^)ii  (at  a  point  close  to  t  =  ^2)  and  of  {k^nh  is  important 
for  the  evolution  of  the  research  program  as  well  as  for  the  assessment  of  the 
validity  of  other  methods  of  data  evaluation.  We  note  in  the  first  place  that  these 
“lower  bound”  and  “true  heat  transfer  coefficients”  are,  respectively,  the  least 
precise  and  least  accurate  values  which  can  be  obtained  from  the  experimental 
data.  Any  conclusion  that  the  precisions  and  accuracies  of  other  methods  of 
data  evaluation  are  lower  than  this  (as  determined  by  repeated  calibrations  for 
“blank  experiments”)  show  that  either  the  experiments  have  been  carried  out 
incorrectly  or  that  the  methods  of  data  evaluation  are  invalid.  We  note,  secondly, 
that  the  Second  Law  of  Thermodynamics  requires  that  >  {k^R)u  (at  the 

same  point  in  time).  Any  conclusion  that  (/c^)ii  >  (^^^2)2  can  only  be  explained 
by  either  errors  in  the  execution  of  the  experiments  or  the  presence  of  variable 
sources  of  excess  enthalpy.  If  the  latter  statement  is  true,  it  is  then  impossible  to 
calibrate  any  conceivable  calorimetric  system  and,  in  that  eventuality,  we  need 
to  rely  on  separate  calibrations  of  the  calorimeters. 


(b)  We  next  carry  out  a  detailed  examination  of  the  (A:^)ii  spreadsheet,  spread¬ 
sheet  6.  We  first  of  all  prepare  a  plot  of  (fei^)ii  ver^  time.  Fig.  40  (compare 
Fig.  4)  on  which  we  also  show  the  derived  values  of  (fc/j)ii.  The  error  bars  show 
±(T  of  (^)ii.  The  “lower  bound  heat  transfer  coefficient”  shows  the  expected 
linear  decrease  with  time  and  the  relevant  regression  line  is  drawn  on  the  plot. 

We  also  examine  the  extrapolations  required  to  obtain  (A:^^)i5i ,  (A:^^)i6i  ?  {k^R  )i7i  ? 
and  Figs.  41  -  44  (compare  Figs.  5-8,  Section  6).  Of  these  plots,  Figs. 

41  and  42  are  probably  soundly  based  in  that  the  origin  of  the  abscissae  (where 
dA9/dt  —  0)  can  be  defined  on  Fig.  38.  The  origin  for  Fig,  43  cannot  be  defined 
and  this  puts  the  mathematical  validity  of  Fig.  44  in  doubt  as  this  requires  a 
combination  of  the  data  used  in  Figs.  42  and  43.  Here,  we  have  said  “prob¬ 
ably  soundly  based”  because  the  question  of  the  origins  for  plots  of  this  kind 
needs  further  investigation  using  appropriate  simulations.  Notwithstanding  our 
reservations  about  the  validity  of  Fig.  44,  we  regard  this  plot  as  the  best  way  of 
estimating  the  “lower  bound  heat  transfer  coefficients”  and  the  water  equivalent 
of  the  cells  in  the  vicinity  of  ^2-  The  water  equivalent  can  also  be  estimated 
from  a  similar  plot  for  {k'^)28i,  see  Fig.  55  below,  and  the  compromise  value  of 
330  JK"^  has  been  used  to  derive  the  (fc^)ii  spreadsheet,  spreadsheet  8.  This 
spreadsheet  is  actually  the  first  iteration  in  the  calculation  scheme.^ ^ 

has  already  been  noted  in  Section  6,  we  do  not  believe  that  this  iteration  has  been 
carried  out  in  the  evaluations  of  the  NHE  data  sets. 


35 


In  view  of  the  scatter  of  the  points  in  Fig.  40,  the  values  of  (fc^^)i6i,  (fc^^)i6ij 
)i7i  (A:^^)i8i  have  not  been  added  to  this  plot  (contrast  the  plot  of 
versus  time,  Fig.  4)  but  are  given  on  a  separate  plot.  Fig.  45,  with  respect  to 
the  regression  line  in  Fig.  40.  Figure  45  also  includes  the  value  of  derived 

by  the  graphical  method,  Fig.  38. 

(c)  In  the  next  step,  we  evaluate  (fc^)i2  and  (fc^)i2  and  plot  these  “true  heat 
transfer  coefficients”  versus  time  in  Fig.  46.  In  carrying  out  this  evaluation, 
we  estimate  the  values  of  Ec  and  6  at  any  given  time  within  ti  <  t  <  t2  hy 
fitting  regression  lines  through  the  data  shown  in  Fig.  38  but  excluding  the 
region  ti  <  t  <  t2  and  those  plots  where  the  temperature  and  cell  potential 
relax  following  perturbations.  It  will  be  evident  t^  the  scatter  of  the  values 
of  (A;^)i2  in  Fig.  46  is  much  larger  than  that  for  (fc^)n  shown  in  Fig.  40  (the 
regression  line  for  (/?^)ii  is  also  shown  in  Fig.  46). 


Figure  47  gives  the  evaluation  for  (A:^^)i62  and  Fig.  48  shows  values  of  (/?^^)i52} 
and  (A:^)2  again  with  respect  to  the  regression  line  for  the 
variation  of  (A:^)ii  with  time,  Fig.  40.  The  increased  scatter  compared  to 
the  data  in  Fig.  45  will  again  be  apparent.  We  conclude  therefore  that  the 
evaluation  of  the  “differential  true  heat  transfer  coefficients”  will  not  be  useful 
for  calculating  the  rates  of  excess  enthalpy  generation.  However,  this  does  not 
mean  that  we  cannot  evaluate  the  differential  rates  -  if  we  should  wish  to  do  so. 
The  way  in  which  we  can  circumvent  the  errors  introduced  by  the  determination 
of  the  “differential  true  heat  transfer  coefficients”  is  discussed  in  subsection  of 
Section  6  (f)  below. 


(d)  We  will  not  describe/discuss  the  use  of  the  first  order  backward  difference  in 
the  evaluation  of  (dAO/dt)  when  considering  the  analysis  of  experimental  data. 


(e)  We  will  consider  next  the  derivation  of  the  various  versions  of  the  integral 
heat  transfer  coefficients.  The  (/c^)2i -spreadsheet,  spreadsheet  7,  gives  the  in¬ 
tegrals  required  for  the  evaluation  of  (A;^)2i  and  (fci^)3i  and  Fig.  49  gives  the 
plots  of  these  coefficients  versus  time.  Figure  50  gives  a  comparison  of  {k'fi)2i 
with  {k'j^)n  and  it  can  be  seen  that  the  use  of  the  integral  procedure  leads  to  a 
marked  reduction  of  the  errors  due  to  the  differentiation  of  “noisy”  experimen- 
tal  data  (which  is  required  for  the  evaluation  of  (fc^)ii  and  the  derived  (A:^)ii). 
Inevitably,  the  values  of  (fc^)2i  and  {k'f^)ii  converge  at  long  times  where  the 
definitions  of  the  integral  and  differential  heat  transfer  coefficients  are  equiva¬ 
lent.  Equally  the  value  of  at  short  times  (as  given  by  the  regression  line 

through  the  points  at  long  times)  converges  onto  the  value  of  {k'j^)ii  at  short 
times  as  these  values  are  again  equivalent,  Fig.  49. 
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It  can  be  seen  from  Figs.  49  and  50  that  the  values  of  (fc^)2i  at  long  times 
deviate  somewhat  from  the  regression  line.  The  reason  for  these  deviations  is 
that  it  is  necessary  to  carry  out  ca  100  integration  steps  in  order  to  suppress 
the  random  errors  in  fi{6)  and  the  enthalpy  input.  Of  course,  the  integration 
procedure  still  leaves  us  with  the  random  errors  in  [9{t)  -  0o]*  However,  it  can 
be  seen  from  Figs.  49  and  50  that  the  further  averaging  to  give  (fc^)2i  has  a 
negligible  effect  on  the  errors  of  this  heat  transfer  coefficient. 

It  can  also  be  seen  from  Fig.  49  that  the  deviations  of  (A;^)3i  at  short  times 
from  the  relevant  regression  line  are  much  larger  than  those  of  at  long 

times.  The  major  reason  for  this  effect  is  that  the  use  of  the  trapezium  rule 
(coupled  to  a  rather  long  measurement  interval)  introduces  appreciable  errors 
into  the  integration  of  /i  (6)  and  the  enthalpy  input  in  regions  where  these  two 
variables  are  changing  rapidly  with  time.  By  contrast,  the  integrations  required 
for  take  place  in  an  initial  region  where  fi{0)  and  the  enthalpy  input 

are  nearly  constant  with  time.  The  curvatures  of  these  variables  with  time  at 
short  times  then  have  negligible  effects  on  the  integrals  and  the  trapezium  rule  is 
perfectly  adequate  for  the  numerical  integrations.  This  was  the  major  reason  for 
our  preference  for  the  use  of  {k'j^)2i  as  compared  to  (fc^)3i-  There  are,  however, 
secondary  reasons  for  this  preference.  These  include  the  fact  that  [6{t)  —  0o] 
is  appreciable  for  the  initial  region  of  the  estimation  of  (^^2)31  whereas  it  is 
negligibly  small  for  the  evaluation  of  {k'ji)2i-  Furthermore,  the  term 

CpM{d-0o)  j  fi{0)dt 

remains  appreciable  throughout  the  range  of  the  evaluation  (so  that  the  values 
of  {k'f^)si  are  subject  to  errors  in  CpM)  whereas  one  can  find  regions  where 

CpMie-do)  I  fiie)dt 

is  zero  (so  that  the  values  of  (fc^)2i  at  these  times  are  independent  of  CpM). 

It  can  be  seen  from  Figs.  49  and  50  that  the  slope  of  the  regression  line  for  the 
variation  of  (A;^)ii  with  time  is  roughly  twice  that  for  the  variations  of  (^^21 
or  (fc^)3i  with  time.  This  is  to  be  expected.  The  time  dependence  of  the  heat 
transfer  coefficients,  Eq.  (21),  does  not  need  to  be  taken  into  account  when 
evaluating  the  differential  version;  it  is  then  the  variation  of  ik'j^)ii  with  time 
which  reveals  this  time  dependence.  This  is  just  the  term  in  Eq.  (21)  and 
we  see  that  we  can  obtain  a  “good”  value  of  7  from  the  slope  of  the  regression 
line,  Fig.  40.  In  the  evaluations  of  (k'j{)2i  and  {k'ji)si  using  Eqs.  (13)  and 
(14)  we  have  usually  defined  fi{6)  using  Eq.  (4)  whereas  our  experience  with 
the  evaluations  of  (A:^)ii  teaches  us  that  we  should  use  the  definition  Eq.  (21). 
Integration  of  this  equation  gives 

{k'n)[j  fi{^)dT --yt  j  fii0)dT +  'y  J  J  fii0)dTdT].  (57) 
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If  we  now  regard  fi{9)  as  being  constant  throughout  the  measurement  cycle 
(which  is  a  rough  approximation  for  the  case  of  the  “lower  bound  heat  transfer 
coefficients”),  then  the  integral  becomes 

(4")/iWi[l-|].  (58) 

It  follows  that  the  heat  transfer  coefficients  given  by  Eqs.  (13)  and  (14)  are 
given  by 

=  +  (59) 

and 

(fcfl)31=(03l[l-Y]  (60) 

within  the  limits  of  this  approximation.  and  (A:^^)3i  are  respectively  the 

values  of  (fc^)2i  at  i  —  T  and  of  (A:^)3i  at  t  =  0.  It  follows  that  the  slopes  of  the 
plots  of  (A;^)2i  and  (A:/i)3i  versus  time  are  one  half  of  the  plot  of  (A;^^)!!  versus 
time. 

Equation  (57)  also  shows  the  way  in  which  we  can  test  whether  the  characteristics 
of  the  Dewar  cells  can  be  described  by  a  single,  time-independent  heat  transfer 
coefficient.  Thus,  evaluation  of  (A:^)2i  according  to  Eq.  (13)  gives  us  the  heat 
transfer  coefficient 

(fck)2i  =  (*fl’)2i  (^)drdr|  /  h  {^)dT  (61) 

so  that  the  “time-independent”  heat  transfer  coefficient  (fc^^)2i  is  readily  deter¬ 
mined.  Figure  51  shows  this  derived  coefficient  versus  the  measurement  interval. 
We  can  see  that  if  we  exclude  the  region  adjacent  to  T  (where  the  methodology 
is  unreliable)  the  values  of  (A;^^)2i  are  within  ±0.01%  of  the  mean  of 
(the  relative  standard  deviation  is  0.0063%).  This  is  the  basis  of  our  statement 
that  the  “integral  lower  bound  heat  transfer  coefficient”  can  be  determined  with 
a  precision  given  by  a  relative  error  of  less  than  0.01%.  We  note  that  the  errors 
in  (A:^^)2i  shown  in  Fig.  51  are  somewhat  larger  than  those  which  we  observed 
in  the  earlier  work. 

The  success  in  deriving  an  unique  value  of  (A:^^)2i  brings  in  its  train  two  further 
aspects.  First  of  all,  we  need  to  assess  the  likely  errors  of  the  “point-by-point” 
values  of  {k'f^)2i  and  (fc^)3i.  Examination  of  spreadsheet  10  shows  that  the 
“minor  term”  ^ 

CpMie-eo)!  [  fiie)dT 

Jt 

is  maximally  ca  ±1.4%  of  the  derived  values  of  (A:^)2i  (at  the  points  t  =  0  and 
i  =  ^2).  On  the  other  hand,  this  “minor  term”  rises  to  ca  25%  of  (A:i^)3i  within 
the  calibration  period  ti  <  t  <  t2  and  is  ca  7%  of  (fc^)3i  at  i  =  0.  Evidently, 
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the  heat  transfer  coefficients  based  on  the  backward  integration  procedure  are 
to  be  preferred  to  those  based  on  forward  integration.  Secondly,  we  see  that 
the  preferred  procedure  would  be  to  use  extrapolation  procedures  to  remove  the 
effects  of  the  water  equivalent,  which  is  the  parameter  subject  to  the  greatest 
uncertainties. 

Figures  52  through  55  show  such  extrapolation  procedures  for  deriving  (A:^'’)25i  > 
(4“)261,  K)271,  and  (fc^^)28i ,  as  well  as  the  corresponding  values  of  CpM. 
The  extrapolation  in  Fig.  52  will  be  to  a  point  in  the  range  0  <  t  <  ti  where 
the  temperature  is  equal  to  that  at  the  end  of  the  measurement  cycle  (roughly 
at  6900s).  As  can  be  seen  the  values  of  and  (A:^^)25i  show  that  this  is 

indeed  so.  The  extrapolations  in  Figs.  53  through  55  should  all  be  to  the  end  of 
the  measurement  cycle  ai  t  =  T.  However,  the  values  of  (A:^^)26i5  (^ij^)27i  and 
(A:^^)28i  are  somewhat  larger  than  that  of  (A:^)2i  at  that  point  (given  the  value 
read  off  from  the  regression  line).  In  1992/93,  we  concluded  that  this  was  due 
to  the  range  of  the  extrapolation  required  being  too  long.  In  order  to  avoid  this 
difficulty,  we  restricted  the  integrations  to  ti  <  t  <t2  and  set  the  origin  of  the 
abscissa  at  a  time  close  to  t  =  ^2-  The  extrapolation,  Fig.  62,  will  now  be  to  the 
value  of  (fcij)ii  at  this  point  and  this  is  again  the  case  (see  also  spreadsheet  8). 

Figure  60  gives  the  “point-by-point”  values  of  (A;J^)2i  for  the  time  range 
<  ^  <  ^2  and  these  are  somewhat  below  the  variation  predicted  based  on 
the  assumption  in  Eq.  (59).  We  note  that  it  has  never  been  resolved  whether 
the  small  discrepancies  observed  are  due  to  using  this  assumption  rather  than 
Eq.  (61)  (however,  they  are  at  least  due  in  part  to  the  use  of  Eq.  (59)  rather 
than  (61)).  Instead,  we  have  recommended  the  use  of  the  extrapolation  proce¬ 
dure  shown  in  Fig.  62.  In  attempting  to  interpret  the  variation  of  {k'^)  with 
time  shown  in  Fig.  60,  it  should  also  be  born  in  mind  that  the  values  of  this 
heat  transfer  coefficient  are  unreliable  for  about  the  first  110  integration  intervals 
adjacent  to  t  =  t2> 

Figure  60  also  gives  the  corresponding  data  for  {k'fi)32  while  Fig.  63  shows  that 
the  extrapolation  procedure  required  for  deriving  (fc^^)36i  is  unreliable,  even  if 
the  first  11  data  points  are  excluded  from  the  extrapolation.  The  main  reason 
for  this  lack  of  reliability  lies  in  the  large  values  of  the  abscissae,  which  reach 
ca  30%  of  the  value  of  (A;^)3i.  (However,  note  that  the  extrapolations  based  on 
the  integrals  spanning  the  whole  of  the  measurement  cycle  0  <  i  <  T,  Figs.  56 
through  59,  are  somewhat  more  satisfactory;  the  procedure  for  deriving  (fc^^)36i  ? 
Fig.  57,  is  reasonably  sound).  It  was  for  this  reason  that  this  procedure  was 
excluded  in  the  Handbook  for  the  ICARUS-1  systems  [1].  We  also  see  that  the 
estimate  of  {k^ii)3i  at  times  close  to  t2  is  somewhat  lower  than  that  predicted 
from  Eq.  (60).  However,  reference  to  spreadsheet  8  shows  that  the  abscissae  are 
still  ca  2.3%  of  the  ordinates  for  the  evaluation  of  {k^ji)3i  at  this  point  (for  the 
24-hour  measurement  cycles  used  in  the  NHE  investigations  this  ratio  rises  to 
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ca  4.4%).  A  3%  error  in  CpM  would  therefore  account  for  such  discrepancies! 


Figure  61  gives  the  values  of  the  “true  heat  transfer  coefficients”  {k‘fi)22  and 
(fc^)32  over  the  restricted  range  h  <  t  <  h  while  Figs.  64  and  65  give  the 
evaluations  of  {k'^)262  and  (A:^^)362*  It  can  be  seen  that  (fc^^)262  is  close  to  the 
value  of  (fei^)ii  at  i  =  ^2  as  given  by  the  relevant  regression  line.  On  the  other 
hand,  (it^^)362  deviates  markedly  from  the  value  of  (fe^)ii  at  t  =  0,  even  if  we 
once  again  exclude  the  first  11  points  from  the  extrapolation.  Spreadsheet  9 
shows  that  the  nominally  “minor  term”  now  reaches  ca  100%  of  the  value  of  the 
’’major  term”  while  the  ’’minor  term”  is  ca  2900%  of  (fe^)32*  It  goes  without 
saying  that  it  is  not  possible  to  evaluate  (fc/e)32  with  this  particular  methodology. 
It  was  for  this  reason  that  the  ICARUS-1  Handbook  recommended  that  {k^ii)32 
should  be  evaluated  “point-by-point”  at  times  close  to  t  —  t2.  Spreadsheet  9 
shows  that  the  “minor  term”  in  the  evaluation  of  {k'j{)32  is  still  ca  12%  of  the 
major  term  while  for  a  24-hour  measurement  cycle  it  would  reach  ca  25%  of 
the  major  term.  The  accurate  determination  of  (^^2)32  is  therefore  fraught  with 
difficulties.  These  conclusions  should  be  compared  with  those  presented  in  the 
poster  given  at  ICCF  7  [4]. 

Finally,  Fig.  66  illustrates  the  determination  of  (A:^^)252j  i-e.,  of  the  “true  heat 
transfer  coefficient”  without  making  use  of  the  calibration  pulse.  Whereas  the  de¬ 
termination  of  the  differential  coefficient  {k^R)i52  fails  (because  of  the  inevitably 
large  errors  introduced  by  the  differentiation  of  “noisy”  experimental  data),  the 
determination  of  (fc^^)252  is  reasonably  successful.  The  evaluation  of  this  version 
of  the  “true  heat  transfer  coefficient”  can  therefore  serve  as  a  useful  check  on 
some  of  the  more  extreme  statements  which  have  been  made  about  the  validity 
of  the  ICARUS- 1  evaluation  procedures. 

(f)  We  come  now  to  the  principal  conclusions  which  we  can  draw  from  the 
detailed  examination  of  “blank”  experiments,  taken  in  conjunction  with  the 
analysis  of  data  generated  by  simulations.  Section  5.  The  optimal  methodology 
for  the  evaluation  of  the  “lower  bound  heat  transfer  coefficient”  is  that  based  on 
the  backward  integration  leading  to  {k'j{)2i,  although  that  based  on  the  forward 
integration  leading  to  (A:^)3i  can  also  be  used.  Furthermore,  the  differential 
form,  (/:^)ii,  is  also  useful. 

However,  the  only  accurate  method  for  the  evaluation  of  the  “true  heat  transfer 
coefficient”  is  that  based  on  the  backward  integrals  and,  especially,  the  extrap¬ 
olation  procedure  giving  (A:^^)262j  which  was  the  methodology  specified  for  the 
ICARUS-1  systems  [1].  It  is  therefore  sensible  to  combine  such  values  of  (fc^^)262 
with  the  corresponding  values  of  (fc^^)26i-  For  the  data  set  discussed  here,  we 
obtain  from  Figs.  60  and  61  (or  Figs.  53  and  64,  spreadsheets  8  and  9): 

-  {k'^)26i  =  0.00043  x  lO-^WR-'*. 


(62) 


This  value  is  to  be  preferred  to 

M^'r)  =  (fck)2  -  =  0.00021  X  10-®WK-S 

given  by  the  graphical  methods  (Eu,  (55)). 


(63) 


It  should  be  noted  that  both  of  these  methods  give  values  of  the  heat  transfer 
coefficients  close  to  the  mid-point  of  the  measurement  cycles  (which  was  one 
reason  for  the  specification  of  the  time  t  =  ^2)-  Furthermore,  the  difference  in 
Eq.  (62)  is  actually  equal  to 

A(fck)  =  ikR)i2  -  ik'R)n  (64) 


at  this  point.  This  is  of  no  particular  consequence  as  far  as  the  evaluation 
of  the  total  excess  enthalpy  for  the  measurement  cycle  is  concerned  because 
the  difference  in  Eq.(62)  applies  to  any  part  of  the  cycle  at  the  first  level  of 
approximation.  It  follows  that 


rO 

Total  excess  enthalpy  =  A (k^)  J  fi(0)dr. 

For  the  particular  example  used  in  the  present  illustration,  we  obtain 


(65) 


Total  excess  enthalpy  ==  102J, 


(66) 


corresponding  to  a  mean  excess  rate  of  0.0006  W. 


However,  we  note  that  the  group  at  NHE  have  attempted  to  calculate  the  vari¬ 
ation  of  the  rate  of  excess  enthalpy  generation  throughout  the  measurement 
cycles.  At  first  sight,  it  would  appear  that  we  need  to  use  the  correct  time- 
dependent  values  of  (A:^)i2,  i*e.,  [^/^(0]i2-  At  the  time  of  writing  of  the  Hand¬ 
book  for  the  ICARUS-1  systems  [1],  it  was  not  clear  how  this  variation  was  to 
be  established.  It  became  clear  subsequently  that  if  the  difference  between  the 
“true”  and  “lower  bound  heat  transfer  coefficients”  could  be  established  at  any 
one  time  (say  A(A:^)t  ),  then  [k'j^{t)]i2  at  any  other  time  would  be  given  by 

[fck(*)]i2  =  [k'Rmhi  +  ^ik'R)M0)tJfiie)t.  (67) 

The  ratio  fi{0)t2  / is  of  order  unity,  which  implies  that  the  correction  term 
is  always  close  to  that  at  the  calibration  point. 


Any  attempt  to  calculate  the  variation  of  rates  of  excess  enthalpy  generation 
within  the  measurement  cycles  must  also  pay  due  regard  to  the  fact  that  it  is 
not  possible  to  calibrate  the  systems  if  the  rate  of  excess  enthalpy  generation 
varies  with  time.  If  that  is  the  case,  then  we  must  derive  A(fc^)  from  separate 
experiments. 
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Equation  (67)  also  points  to  a  further  important  conclusion.  Again  at  the  time 
of  writing  of  the  Handbook  for  the  ICARUS-1  systems  [1],  we  believed  that  the 
precision  of  (k'fi)  12  (and  of  other  versions  of  the  “true  heat  transfer  coefficient”) 
would  always  be  given  by  the  accuracy  of  that  coefficient  which  is  certainly  lower 
than  the  precision  of  (/i:^)ii.  Equation  (67)  shows  that  this  is  incorrect.  The 
precision  of  (fc/^)i2  is  nearly  identical  to  the  precision  of  (fei^)ii.  It  follows  that 
changes  in  the  rates  of  excess  enthalpy  generation  can  be  established  at  the  same 
level  of  precision  as  that  of  i.e.,  with  relative  errors  ca  0.01%. 
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SECTION  7:  ASSESSMENT  OF  THE  SPECIFICATION  OF  THE 
ICARUS-1  EXPERIMENTAL  PROTOCOLS  AND  DATA  EVALU¬ 
ATION  PROCEDURES. 

The  specification  of  the  ICARUS-1  experimental  protocols  and  data  evaluation 
procedures  has  been  outlined  in  Section  4,  It  is  now  important  to  assess  the 
usefulness  and  validity  of  this  specification  in  the  light  of  its  application  to  data 
generated  by  simulations  (see  Section  5)  and  by  “blank”  experiments  (see  Section 
6).  Of  course,  this  procedure  is  “back-to-front” :  the  specification,  Section  4,  was 
evolved  from  a  consideration  of  the  type  of  results  outlined  in  Sections  5  and  6. 


It  is  also  important  to  consider  the  “raison-d’etre”  of  this  part  of  the  research 
program  and  some  comments  on  this  are  given  at  the  end  of  this  section.  Fur¬ 
thermore,  it  is  important  to  consider  the  causes  of  failures  to  achieve  evaluations 
especially  in  the  light  of  the  well  publicized  publication  from  the  group  at  NHE 
[5].  Finally,  it  is  necessary  to  consider  “short-cuts”  to  achieving  satisfactory 
data  evaluations. 

In  our  view,  the  major  results  derived  from  Section  6  (backed  up  by  the  in¬ 
vestigation  in  Section  5)  are  the  comparisons  of  and  (fc^^)!!  in  Fig.  50 

and  the  reduction  of  (A:^)2i  to  a  single,  time-independent,  lower  bound  heat 
transfer  coefficient.  Fig.  51.  The  comparison.  Fig.  50,  illustrates  immediately 
the  need  to  avoid  the  differentiation  of  “noisy”  experimental  data  (required  for 
the  evaluation  of  (A:^)ii  and  the  benefits  of  using  instead  the  integration  pro¬ 
cedures  in  deriving  (A:^)2i);  however,  see  further  below.  This  was  the  basis  for 
the  specification  of  the  construction  of  the  ik'j^)2i  spreadsheets  following  on  the 
construction  of  the  (fc^)ii  spreadsheets  (see  (vi)  and  (vii)  of  Section  4). 

Results  such  as  those  illustrated  in  Fig.  50  show  that  it  is  possible  to  interpret 
the  systematic  variations  with  time  of  ca  0.4%  of  the  “integral  lower  bound  heat 
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transfer  coefficient”  while  Fig.  51  shows  that  it  is  possible  to  reduce  such  data  to 
a  single,  time-independent,  heat  transfer  coefficient,  (fc^^)2i  with  relative  errors 
below  0,01%.  This  result  is  hardly  surprising.  The  “physics”  of  the  calorimeters 
are  quite  simple  (they  are  “ideal  well  stirred  tanks”)  and  the  errors  are  mainly 
due  to  those  set  by  the  temperature  measurements.  It  is  also  relatively  straight¬ 
forward  to  specify  the  changes  which  would  need  to  be  made  to  reduce  the  errors 
“  say,  to  0.001%  -  if  that  should  ever  prove  to  be  necessary  or  desirable. 

The  comparison  of  (fc^)2i  with  (fc^)3i  in  Fig.  50  as  well  as  in  spreadsheet  7 
shows  that  the  errors  in  (A;/^)3i  are  inevitably  larger  than  those  of  This 

is  mainly  due  to  the  larger  contribution  of  the  term 

CpM{e-eo)/  f  fimdr 

Jo 

to  (/c^)3i  rather  than  the  corresponding  contribution  of 

CpMie-do)/  J^fii0)dT 

to  (A:^)2i-  We  can  see  immediately,  that  given  the  option  of  using  forward 
or  backward  integration  we  should  use  the  latter  in  order  to  achieve  accurate 
evaluations.  This  restriction  becomes  much  more  important,  however,  when  we 
consider  the  derivation  of  the  “true  heat  transfer  coefficients,”  (fci^)22  and  (A:^)32  ? 
as  well  as  of  the  likely  effects  of  timing  errors  on  the  evaluations.  Such  evaluations 
must  necessarily  be  mainly  restricted  to  the  duration  of  the  calibration  pulse, 
<  t  <  t2^  Figures  60  and  61  and  the  associated  spreadsheets  8  and  9  show 
that  whereas  an  initial  assessment  might  well  be  based  on  the  use  of  (A:^)3i5  the 
evaluation  and  use  of  the  “true  heat  transfer  coefficient,”  (/c^)32,  must  follow 
strictly  the  instructions  laid  down  in  the  ICARUS-1  Handbook  [1]  (see  (vii)  of 
Section  4)  -  and,  even  then,  the  errors  are  much  larger  than  those  of  {k'j^)22- 

The  optimum  methodology  for  deriving  the  diiference  between  the  “true”  and 
“lower  bound  heat  transfer  coefficients,”  i.e.,  A(A:^),  Eq.  (62),  is  to  estimate 
(4^)262  from  plots  such  as  that  in  Fig.  64  and  {k'^)26i  from  the  correspond¬ 
ing  plot  in  Fig.  62.  Such  evaluations  automatically  eliminate  the  contributions 
of  the  water  equivalent  on  the  estimates,  a  parameter  which  is  subject  to  the 
greatest  degree  of  uncertainty.  The  evaluation  of  the  difference  A(A:^)  allows 
us  also  to  circumvent  the  errors  in  the  evaluation  of  the  “true  differential  heat 
transfer  coefficient,”  ik'f{)i2i  Fig.  46,  because  A(/;:^)  applies  equally  to  (A:^)i2 
and  {k'f^)u,  Eq.  (64),  as  to  (fc^^)262  and  (fc^^)26ij  Eq,  (62).  It  follows,  therefore 
that  we  can  evaluate  rates  of  excess  enthalpy  generation  at  the  level  of  accuracy 
determined  by  A(/c^)  and  precision  of  (A:^)ii.  However,  we  cannot  see  why  we 
should  ever  wish  to  do  so,  given  that  we  have  had  to  evaluate  the  backward 
integrals  to  determine  A(fc^)  so  that  the  values  of  the  total  excess  enthalpy  in 


a  measurement  cycle  (or  any  given  part  of  a  measurement  cycle)  follow  imme¬ 
diately  as  by  the  application  of  Eq.  (65).  If  we  evaluate  instead  the  rates  of 
excess  enthalpy  generation,  then  we  must  sum  these  rates  (multiplied  by  the 
measurement  interval)  to  obtain  the  total  excess  enthalpy  [4]  (more  exactly,  we 
need  to  apply  an  appropriate  integration  rule). 


transfer 


R  )l62, 


Sections  5  and  6  also  illustrate  the  evaluations  of  the  ’’time-independent”  heat 
coefficients  (fc^^)i5i,  (A:^^)i6i,  (fc^^)i7i,  (*^^)i8i,  {k^R)ih2, 

(fc^^)2ol?  (fc^^)261^  (fcj4^)271j  (fcif)281>  (fc^^)252  5  {k^  )262j 
(fefi^)272,  (fc^^)3ol,  (^ii°)361?  (fei?^)37l5  {k'j^)381^  (fc^^)352,  (fcfi^)362  and 

{k^^)s72-  Of  these,  the  evaluations  of  (fc^^)26i  and  (/c^)262,  Figs.  62  and  64, 
were  the  methodologies  specified  for  the  ICARUS-1  system.  The  reduced  preci¬ 
sions  and  accuracies  of  the  series  {k^^)ij,i  compared  to  those  of  {k'^)2j,i  will  be 
apparent  as  will  be  the  reduced  precisions  and  accuracies  of  the  series  {k'j{)3,j.i 


compared  to  those  of  {k'^)2jj-  The  evaluation  of  {k'^)sQ2y  Fig.  65,  is  especially 
prone  to  error.^^ 


The  failure  to  achieve  satisfactory  calibrations  of  the  cells  using  the  procedures 
leading  to  (A:^)362  in  the  studies  by  the  group  at  NHE  coupled  to  the  evident 
timing  errors  in  the  ICARUS-1  system  installed  in  Sapporo,  also  prompted  an 
investigation  of  the  likely  effects  of  such  errors  if  these  were  not  properly  taken 
into  account.  Figs.  26-37  show  the  expected  results.  Whereas  the  “target 
evaluations”  of  {k'^)26i  and  (fc^°)262  are  relatively  insensitive  to  the  effects  of 
such  errors  (Figs.  26-29),  the  evaluations  of  (fc^°)36i  and  (fc^^)362  are  markedly 
degraded  (Figs.  30-37).  Moreover,  the  values  of  {k'^)3Gi  and  {k^^)3G2  (the  inter¬ 
cepts  of  the  plots)  now  depend  on  the  time  range  of  the  regression  fit  (the  dis¬ 
cussion  of  these  effects  is  beyond  the  scope  of  the  present  report).  It  is  evidently 
critically  important  to  avoid  the  effects  of  such  timing  errors  because  there  are 
other  reasons  which  preclude  the  achievement  of  satisfactory  calibrations  most 
notably  the  effects  of  changes  in  the  rates  of  excess  enthalpy  generation  espe¬ 
cially  those  due  to  “positive  feedback,”  e.g.,  see  [4,12].  It  is  essential  therefore 
that  tests  of  the  performance  of  the  instrumentation  be  carried  out  using  “blank 
experiments”  (see  Section  4  (iii))  in  order  to  avoid  the  complications  introduced 
by  the  study  of  the  Pd-D20  system.^^ 


Sections  5  and  6  will  also  have  illustrated  the  reasons  for  the  need  to  use  48- 
hour  measurement  cycles  [see  Section  4  (i)]  as  well  as  the  benefits  of  examining 
and  evaluating  plots  of  the  raw  data  such  as  that  in  Fig.  38  [see  Section  4  (v)]. 

^"^Unfortunately,  such  evaluations  appear  to  have  been  a  major  part  of  the  investigations 
carried  out  by  NHE. 

^®We  have  never  been  able  to  obtain  data  for  the  study  of  any  such  “blank  systems”  which 
may  have  been  carried  out  by  NHE. 
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Finally,  it  is  important  to  assess  the  usefulness  of  the  calibration  procedures 
described  in  this  section  as  well  as  the  preceding  sections.  It  appears  that  this  is 
a  topic  which  has  been  misunderstood.  It  would  hardly  be  possible  to  investigate 
all  measurement  cycles  of  all  experiments  at  the  level  of  detail  set  out  in  Section 
6.  The  real  purpose  of  developing  precise  and  accurate  methods  of  calibration 
is  to  gain  an  adequate  understanding  of  the  “instrument  function”  and  then  to 
use  the  calibration (s)  to  measure  the  rates  of  excess  enthalpy  generation,  A 
secondary  objective  (related  to  the  assessment  of  the  calibrations)  is  to  point 
the  way  towards  necessary  (or  useful)  improvements  of  the  instrumentation.  We 
note  here  that  the  original  cells,  which  were  not  silvered  in  the  top  section,  gave 
a  flat  temperature-time  base  line  because  changes  in  the  heat  transfer  coefficient 
with  time  were  compensated  by  changed  in  the  enthalpy  input.  However,  the 
heat  transfer  coefficients  varied  markedly  with  time  and  this  complicated  the 
methods  of  data  analysis.  The  silvering  in  the  top  sections  of  the  ICARUS- 
1  and  -2  calorimeters  has  markedly  reduced  this  time  dependence  of  the  heat 
transfer  coefficients,  but  the  effects  (e.g.,  see  Fig.  50)  do  still  have  to  be  taken 
into  account.^®  In  1990/91,  it  appeared  important  that  we  should  investigate 
a  further  modification  of  the  calorimeter  design  designated  as  the  ICARUS-4 
version  in  1993  (the  ICARUS-3  designs  were  never  constructed)  but  redesignated 
since  that  time  as  ICARUS-14,  vol.  I,  Fig.  27.  We  believe  that  this  design  would 
give  a  flat  base  line  for  the  heat  transfer  coefficient-time  plot  and  that  this  would 
also  make  the  heat  transfer  coefficient  insensitive  to  the  operating  conditions. 


An  alternative  approach  would  be  to  investigate  the  maintenance  of  constant 
electrolyte  levels  as  has  been  done  in  the  investigations  carried  out  by  the  group 
at  Grenobles  [15].  It  would  be  important  to  try  to  secure  the  release  of  data  for 
“blank  experiments”  carried  out  by  that  group. 


A  secondary  reason  for  carrying  out  repeated  calibrations  is  to  monitor  the 
system  behavior  [4,  5]  especially  the  effects  of  the  onset  of  positive  feedback. 
However,  such  investigations  do  not  need  to  be  carried  out  at  the  high  levels  of 
precision  and  accuracy  required  for  the  calibration  of  the  calorimeters. 

note  that  this  has  not  been  taken  into  account  in  the  work  carried  out  by  the  group 
at  NTIE. 
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Fig.  2.  ICARUS-2  simulation.  Temperature-time  and  input  enthalpy-time  series. 
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Fig.  4.  ICARUS-2  simulation,  (k'c)^,  (k'c°)i5i,  (k'c°)i6i.  (k'c°)i7i.  and  (k'c°)i8i  versus  time. 
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Fig.  6.  ICARUS-2  simulation.  Evaluation  of  (k'c°)i6i 
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Fig.  8.  ICARUS-2  simulation.  Evaluation  of  (k'c°)i8i 
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Fig.  10.  ICARUS-2  simulation.  Evaluation  of  (k'c°)i62- 
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Fig.  12.  ICARUS-2  simulation.  Evaluation  of  (k'c°)i5i  using  first  order  differences. 
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Fig.  16.  ICARUS-2  simulation.  (k'c)i2,  {k'c)i52.  (•<'c°)i62-  ^i^cl  (k'c)i72  versus  time  first  order  differences. 
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Fig.  18.  ICARUS-2  simulation.  (k'c)2i  and  (k'c)3i  versus  time. 
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Fig.  22.  ICARUS-2  simulation.  Evaluation  of  (k'cO)36i 


Fig.  24.  ICARUS-2  simulation.  Evaluation  of  (k'c°)36i  using  mid-point  integration  rule. 


Fig.  26.  ICARUS-2  simulation.  Evaluation  of  (k'RO)26i  with  +30J  error  in  input  enthalpy. 
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Fig.  28.  ICARUS-2  simulation.  Evaluation  of  (k'pjO)26i  with  +30J  error  in  input  enthalpy. 
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Fig.  30.  lCARUS-2  simulation.  Evaluation  of  (k'c°)36i  +30J  error  in  input  enthalpy. 
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Fig.  32.  ICARUS-2  simulation.  Evaluation  of  (k'c°)362  with  error  +30J  in  input  enthalpy  and  using  intervals  156-177. 
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Fig.  34.  ICARUS-2  simulation.  Evaluation  of  (k'cO)36i  with  error  +30J  in  input  enthalpy. 


Fig.  36.  ICARUS-2  simulation.  Evaluation  of  (k'cO)362  with  error +30J  in  input  enthalpy  and  using  intervals  156-177. 
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Fig.  38.  The  "raw  data"  for  a  blank  experiment  (Pt  in  0.1  M  Li0D/D20) 
in  an  ICARUS-2  calorimeter.  Third  measurement  cycle.  Evaluation  of 
(k'4)Hi  and  (k'D)2  using  the  graphical  method. 
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Fig.  40.  (k'R)ii  for  the  experiment  iliustrated  in  Fig.  38  plotted  versus  time. 
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Fig.  44.  ICARUS-2.  Evaluation  of  (I<'r°)i8i- 
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Fig.  56.  Evaluation  of  (k'R°)35i. 
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Fig.  62.  ICARUS-2.  Evaluation  of  (k'R°)26i  using  the  ICARUS-1  Method 
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Fig.  63.  ICARUS-2.  Evaluation  of  (k'R°)36i  using  the  ICARUS-1  Method. 
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Fig.  64.  Evaluation  of  (k'FiO)262  using  the  ICARUS-1  Method 
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Fig.  66.  Evaluation  of  (k'R°)252  using  the  ICARUS-1  Method. 
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